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ABSTRACT 

An  experiment  was  performed  with  broad  spectrum  acoustic  waves  . 
Continuous  broad-spectrum  illumination  was  employed  to  examine  the  fre¬ 
quency  dependence  of  backscatter  from  statistically  rough  surfaces.  In 
addition,  the  operation  of  a  panchromatic  radar  system  was  simulated. 

The  frequency  dependence  of  backscatter  was  measured  for  an  undu¬ 
lating  surface  appropriate  for  application  of  the  Kirchhoff-Huygen's  Principle, 
and  for  a  very  rough  surface.  The  frequency  dependence  of  the  very  rough 
surface  was  found  to  be  in  agreement  with  experimental  results  of  other 
programs.  The  frequency  dependence  of  the  undulating  surface  was  compared 
with  the  theory  based  on  the  Kirchhoff  Method  and  was  found  to  agree  only 
in  certain  frequency  regions.  The  results  support  the  contention  that  the 
effective  surface  statistics  are  both  frequency  and  angle  dependent  and  that 
the  dependences  are  interrelated.  However,  theoretical  treatment  was 
found  that  incorporated  the  correct  dependences  to  predict  the  frequency 
dependence  measured  in  this  experiment. 

In  recent  years  the  use  of  imaging  radar  systems  for  both  military 
reconnaissance  and  geoscience  investigations  has  expanded  rapidly. 
Associated  with  this  expanding  demand  has  been  the  growing  realization 
that  the  monochromatic  nature  of  such  system  handicaps  their  utility, 
especially  in  military  reconnaissance  applications.  This  work  is  an  experi¬ 
ment  with  a  panchromatic  system. 

A  system  has  been  designed  and  constructed  to  implement  a  pan¬ 
chromatic  technique  using  acoustic  waves  in  water.  The  validity  of  using 
acoustic  simulation  of  electromagnetic  waves  has  been  well  established, 
and  the  advantages  of  doing  so  are  clearly  evident  in  this  study.  An  evalua- 


tion  of  the  panchromatic  system  and  verification  of  its  utility  are  presented. 

Images  were  obtained  of  object  targets  using  monochromatic  and 
swept-frequency  illumination.  The  images  clearly  illustrate  the  effects 
of  scintillation  and/or  glint  on  monochromatic  images,  and  the  elimination 
of  the  effect  gained  by  using  panchromatic  illumination.  The  techniques 
proved  extremely  effective  for  averaging  the  lobes  of  complex  reradiation 
patterns  and,  hence  improving  the  target  definition .  This  improvement 
was  gained  by  employing  less  than  +  6%  frequency  bandwidth.  The  results 
indicate  than  panchromatic  systems  should  offer  significant  advantages  over 
monochromatic  systems  in  applications  such  as  terrain  imaging,  altimetry, 
and  tracking  systems. 
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CHAPTER  1 


INTRODUCTION 


An  experiment  was  performed  with  bread- spectrum  acoustic  waves. 
Continuous  broad- spectrum  illumination  was  employed  to  examine  the 
frequency  dependence  of  backscatter  from  statistically  rough  surfaces.  In 
addition,  the  operation  of  a  panchromatic  radar  system  was  simulated. 

Theoretical  investigations  of  electromagnetic  and  acoustic  wave 
scattering  have  been  conducted  for  many  years  (Beckmann  and  Spizzichino, 
1963;  Fung,  1966;  Daniels,  1963;  Parkins,  1966;  Rayleigh,  1896;  Rice,  1951; 
Bass,  1961;  Ament,  et  al.  1959).  Unfortunately  the  extent  of  experimental 
verification  available  for  carefully  measured  surfaces  is  extremely  limited 
(Earing  and  Smith,  1966).  In  this  work  comparison  has  been  made  between 
experimental  and  theoretical  descriptions  of  the  frequency  dependence  of 
backscatter  from  randomly  rough  surfaces. 

Radar  has  been  employed  for  many  years  in  detection  and  ranging 
applications  (Skolnik,  1962;  Barton,  1964;  Kerr,  1951).  More  recently, 
attempts  have  been  made  to  utilize  radar  in  geoscience  applications  that 
require  much  more  information  from  the  return  signal  than  the  existence  and 
range  of  a  target  (Simonett,  1966;  Pierson,  etal.,  1965;  Moore,  1962; 
MacDonald,  et  al.  1968).  Much  of  the  desired  information  is  believed  to 
be  obtainable  by  observing  the  angle,  polarization,  and  frequency  dependence 
of  the  illuminated  surface  (Rouse,  etal.,  1966;  Ellermeier,  etal.,  1966; 
Peake,  1958).  Some  experimental  radar  systems  have  been  used  to  measure 
angular  dependence  of  backscatter  (Ament,  et  al.,  1959;  Ohio  State  Univer¬ 
sity,  1963).  Likewise  some  limited  applications  of  the  polarization  depen¬ 
dence  have  been  employed  (Moore  and  Dellwig,  1967;  Katz  and  Spetner,  1960) 
Polychromatic  (multi- frequency)  systems  are  also  not  new  (Grant  and  Yaplee, 
1957).  Unfortunately,  most  experimental  results  presently  available  are  in¬ 
adequate  for  use  in  verifying  the  theoretical  predictions  of  these  dependencies 
Polychromatic  systems  have  been  used  for  measurements  of  back¬ 
scatter  from  object  targets  and  terrain  (Ohio  State  University,  1963;  Ament, 
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etal.,  1959;  Grant  and  Yaplee,  1957;  Wiltse,  etal.,  1957).  However, 
their  fundamental  purposes  have  been  to  provide  data  for  either  determination 
of  an  appropriate  design  irequency  for  monochromatic  radar  systems  used  for 
particular  applications  or  design  information  needed  to  determine  necessary 
loop  gain  for  particular  monochromatic  radar  systems.  The  added  information 
about  the  target  characteristics  intuitively  believed  to  be  present  in  the 
behavior  of  the  return  at  different  frequencies  has  not  been  extracted  except 

to  an  extremely  limited  extent  (Katz,  1966). 

Radar  systems  incorporating  some  degree  of  frequency  agility  have 
been  constructed  for  applications  requiring  anti- jamming  capability  or  a 
reduction  of  the  effects  of  fading  or  scintillation  (Skolnik,  1962).  These 
systems  are  not  appropriately  described  as  panchromatic  since  the  frequency 
spread  is  usually  very  small,  however,  as  is  discussed  in  Chapter  4  and  6, 
the  primary  utility  of  panchromatic  systems  is  the  reduction  of  the  effect  of 
positional  fading  or  scintillation.  In  general  the  frequency  dependence  of 
backscatter,  at  least  from  terrain  targets,  is  sufficiently  broad  that  the  pan¬ 
chromatic  system  does  not  acquire  frequency  "signature"  information  as  such 

(Katz,  1966). 

Polychromatic  and  polypanchroma  tic  (multi- panchromatic)  systems 
should  permit  determining  frequency  "signatures"  for  object  targets  and 
terrain.  The  discussion  of  polychromatic  radar  experiments  in  Chapter  2 
and  the  results  of  the  acoustic  experiment  presented  in  Chapter  5  indicate 

the  information  expected  from  such  systems. 

The  experiments  were  conducted  using  acoustic  waves  in  water. 

It  is  important  that  the  reader  appreciate  the  analogy  between  acoustic 
waves  and  electromagnetic  waves,  and  understands  that  the  results  obtained 
by  this  research  are  applicable  to  broad- spectrum  radar  systems.  Specific 
cases  where  the  electromagnetic  wave  and  acoustic  wave  analogy  is  imper¬ 
fect,  such  as  scattering  involving  polarization  effects,  are  pointed  out  in 
the  text;  however  in  general  one  can  assume  that  the  results  obtained  by 
measurements  of  acoustic  wave  returns  from  a  target  may  just  as  well  have 
been  obtained  by  flying  a  radar  over  a  particular  terrain  region  having  the 

same  statistical  properties. 

Acoustic  simulation  of  electromagnetic  problems  is  well  developed 
and  documented.  Even  so  its  value  as  an  analog  computer  is  either  not 
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understood  or  not  appreciated  by  many  radar  engineers.  An  acoustic  simu¬ 
lator  makes  it  possible  to  perform  in  the  laboratory  many  of  the  experiments 
normally  reserved  for  flight  testing,  and  to  adjust  parameters  of  the  radar 
system  in  a  highly  controlled  and  repeatable  manner .  The  simulator  is 
advantageous  because  both  dista  ce  and  time  are  scaled  by  a  factor  contain¬ 
ing  the  ratio  of  acoustic  propagation  in  water  to  the  electromagnetic  propaga¬ 
tion  in  air;  a  factor  of  2  x  105.  Therefore,  the  radar  problem  is  compressed 
to  a  size  easily  accommodated  in  a  water  tank,  and  the  electronic  system 
operates  in  the  region  near  1  MHz,  rather  than  at  microwavelengths.  The 
flexibility  is  further  enhanced  by  the  option  to  scale  distance  in  a  non-linear 
manner.  Hence,  model  experiments  can  be  performed  rapidly,  repeatably,  and 
at  much  lower  cost  than  full  scale  radar  experiments,  and  can  provide  results 
completely  analogous  to  a  large  number  of  electromagnetic  problems  (Edison, 
1961).  A  detailed  discussion  of  the  analogy  and  its  validity  is  contained  in 
Chapter  4. 

The  design  and  construction  of  a  system  providing  broad- spectrum 
illumination  occupied  much  of  the  effort  in  this  program.  Of  particular  con¬ 
cern  is  the  behavior  of  the  ultrasonic  transducers  used  as  the  transmitter 
and  receiver  antennas.  These  devices  are  piezoelectric  crystals  which  are 
resonant  at  particular  discrete  frequencies.  Operation  of  these  devices 
off- resonance,  as  is  necessary  in  the  broad- spectrum  mode,  poses  special 
problems  because  of  the  sharp  skirts  of  the  transducer  response  and  spurious 
resonances  due  to  the  transducer  mountings.  Overcoming  these  problems 
and  defining  the  transducer  effects  on  the  data  are  significant  achievements 
of  the  system  effort.  A  thorough  discussion  of  the  transmission  and  recep¬ 
tion  methods,  data  processing,  and  calibration  is  necessary  for  interpretation 
of  the  experimental  results  and  establishing  their  validity.  This  discussion 
is  Chapter  3 . 

This  research  program  is  divided  into  two  sections:  1)  scattering 
experiment  and  2)  image  experiment.  The  scattering  experiment  consisted 
of  the  measurement  of  backscattering  from  randomly  rough  surfaces;  the  deter¬ 
mination  of  the  frequency  dependence  predicted  by  scattering  theory;  and  the 
comparison  of  experimental  and  theoretical  results.  Measurements  of  back- 
scatter  power  from  randomly  rough  surfaces  were  recorded  for  incidence  angles 
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from  vertical  to  60  degrees  across  a  frequency  range  of  about  30:1  (\=  15mm 
to  0.5mm  in  water).  A  scattering  surface  was  used  that  is  gently  undulating, 
with  slopes  large  relative  to  the  incident  wavelengths.  In  addition  measure¬ 
ments  were  made  using  the  basic  surface  with  small  scatterers  superimposed. 
The  scattering  experiment  procedure  and  results  are  presented  in  Chapters 
4  and  5  respectively. 

The  image  experiment  consisted  of  recording  in  image  form  (B-scan) 
the  return  from  particular  object  and  area  extensive  targets.  The  principal 
objective  was  to  examine  the  effect  of  broad- spectrum  illumination  on  the 
return.  The  monochromatic  reradiation  lobing  structure  of  the  targets  was 
recorded  and  the  "blending"  characteristic  of  the  broad- spectrum  illumination 
was  introduced  and  recorded.  The  image  experiment  results  are  presented 
in  Chapter  6 . 
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CHAPTER  2 


BACKGROUND 


The  study  of  waves  scattered  from  rough  surfaces  has  been  active 
for  several  decades  because  of  the  need  of  solutions  to  application  prob¬ 
lems .  The  problems  of  primary  concern  have  been  associated  with  radar, 
but  the  studies  relate  to  line-of-sight  radio  communication,  radar  astro¬ 
nomy,  and  sonar  among  others  (Beckmann  and  Spizzichino  1963).  Although 
a  general  and  exact  solution  to  this  problem  is  yet  lacking ,  considerable 
progress  has  been  made  in  developing  various  special  methods  that  are 
valid  in  many  situations. 

In  recent  years  radar  has  found  utility  as  a  sensor  for  geoscience 
investigations  (Pierson,  et  al.  1965).  Such  applications  require  that  the 
radar  return  provide  information  concerning  the  nature  of  the  illuminated 
surface;  for  example,  its  topography,  moisture  content,  surface  cover, 
etc.  The  radar  parameters  of  use  in  describing  the  desired  terrain  proper¬ 
ties  are  incidence  angle,  polarization,  and  frequency  (Rouse,  et  al.  1966). 
Therefore,  at  least  for  radar  problems,  expressions  derived  for  average 
backscatter  power  must  correctly  incorporate  the  dependence  of  the  ampli¬ 
tude  on  these  parameters. 

Many  authors  have  treated  the  rough  surface  scattering  problems 
by  statistical  methods  (Fung,  1966).  Expressions  have  been  derived  which 
attempt  to  incorporate  the  angle,  polarization,  and  frequency  dependence 
of  the  mean  return  power  (Ohio  State  University,  1963).  The  Kirchhoff  method 
is  especially  popular  for  this  effort.  The  expressions  obtained  by  this 
approach  incorporate  terms  representing  the  statistics  of  the  surface. 

Attempts  have  been  made  to  obtain  the  statistical  parameters  of  a  surface,  the 
moon  surface  for  example,  by  fitting  the  equations  to  the  backscatter  data 
at  difierent  incidence  angles  and  wavelengths  (Evans,  1957;  Hayre ,  1964). 

It  has  been  found  that  the  values  for  the  correlation  distance  and  the  standard 
deviation  dre  different  for  different  wavelength  incident  signals  (Fung,  1965). 
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The  conclusion  has  been  that  such  parameters  are  wavelength  dependent 
and  that  the  incident  wavelength  acts  as  a  "size  filter". 

2.1  Theoretical  Analyses 

A  direct  attempt  to  explain  the  wavelength  dependence  of  scatter 
was  made  by  Spetner  and  Katz  (i960)  using  two  different  scattering  models. 
Tneir  results  are  presented  below.  Chapter  5  offers  another  explanation 
employing  expressions  obtained  by  the  Kirchhoff  method,  which  is  summar¬ 
ized  below  and  detailed  in  Appendix  A. 


2.1.1  Soetner  and  Katz  Random  Scatterer  Model 


The  surface  being  considered  is  assumed  to  be  composed  of  a  large 
number  of  incoherent  independent  scatters.  The  normalized  radar  cross 
section  of  the  surface  is  the  product  of  the  density  of  scatterers  per  unit 
area  and  the  average  cross  section  of  a  single  scatterer.  For  a  sufficiently 
short  incident  wavelength  the  individual  scatterers  are  distinct.  However, 
as  \  increases  adjacent  scatterers  tend  to  merge  and  appear  as  a  single 
scatterer,  hence  the  effective  density  of  scatterers  decreases.  This  is 
expressed  as 


X  <  X, 
X  >  X, 


(2.1) 


where  p^  -  density  of  scatterers  per  unit  area 

A,  =  wavelength  above  which  scatterers  become  coherent 

A  =  incident  wavelength 

n  C  -  constants 
ro>  1 

The  average  scatterer  cross  section  for  small  wavelengths  (but  not 
necessarily  as  small  as  ) ,  is 


(2.2) 
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where 


/o  =  power  reflection  coefficient 
/4  =  effective  area  of  the  scatterer 


For  large  wavelengths  a  Rayleigh-type  scattering  is  assumed  and  the 
scattering  cross  section  is 


where  V  =  volume  of  the  scatterer 

C ^  =  constant 

The  transition  wavelength  between  the  application  of  (2.2)  and  (2.3)  is 
denoted  Xq  .  The  results  of  this  approach  are  summarized  in  figure  2.1. 


2.1.2  Spetner  and  Katz  Specular  Point  Model 


The  surface  being  considered  is  reflecting,  is  of  irregular  shape 
but  continuous,  and  has  continuous  derivatives.  The  contributions  to  the 
radar  return  are  assumed  to  come  predominantly  from  two  sources:  the 
energy  which  is  reflected  from  large  facets  of  the  surface  which  are  oriented 
perpendicular  to  the  line-of-sight,  and  the  energy  which  is  scattered  iso¬ 
tropically  by  the  small  scatterers  . 

The  scattering  coefficient  is  assumed  to  be  of  the  form 


(2.4) 


where 


r-%  =  scatterer  densities  for  small  and  large  scatterers 

respectively . 

O’  cr  -  scattering  cross  sections  for  small  and  large 
scatterers  respectively. 


The  dependence  of  Pg  on  X  is  the  same  as  shown  in  equation  (2.1). 
The  function  p^.  expresses  the  density  of  facets  oriented  normal  to  the 
incident  energy.  Therefore  it  is  related  to  the  slope  spectrum  of  the  surface; 


(DB)  (arbitrary  zero)  <j°  (DB)  (arbitrary  zero) 


2 

in  particular  it  is  related  to  the  slope  variance,  s  .  Assuming  a  flat-slope 
spectrum  with  a  cutoff  at  X9  'L ,  and  considering  the  slope  variance  to  be  an 
effective  parameter  which  is  related  to  the  measured  variance,  Sq,  by 


then  the  density  of  large  scatterers,  at  vertical  incidence,  is 


(2.5) 


The  scattering  cross  sections  are 


X  <  A* 
A  >  X* 


(2.6) 


and 


(flat  plate) 


(2.7) 


(Rayleigh  scatter)  (2.8) 


The  final  result  obtained  for  a  surface  composed  as  described  and  having 
a  Gaussian  slope  distribution  and  flat-slope  spectrum  is 


C-coi  Ze/z  s?) 


;  A  <  A* 

(2.9) 


*  S 

+  j  A  > 
A4 


The  results  for  C5  <  C3  C4  are  shown  in  figure  2.2. 
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The  concepts  advanced  in  this  work  are  basically  that  various  sub¬ 
ranges  of  scatter  sizes  effect  the  return  differently  at  different  wavelengths, 
thus  the  average  return  power  expression  must  somehow  accept  the  role  of 
the  changing  range  of  dominant  scatterers  as  the  wavelength  varies.  This 
idea  is  elaborated  in  section  5.4  except  that  the  Kirchhoff  method  is 
employed . 

2.1.3  Kirchhoff  Method 

The  method  of  physical  optics,  or  Kirchhoff  method,  is  based  on 
the  assumption  that  the  incident  field  is  "loca^xy  "  reflected  at  a  point  on 
the  irregular  surface  as  if  an  infinite  plane  were  tangent  to  the  surface  at 
that  point.  This  assumption  does  not  restrict  the  method  to  roug'.ness  of 
any  particular  size  relative  to  the  wavelength.  However,  its  use  is  restricted 
to  surfaces  which  are  "locally  flat"  with  respect  to  wavelength.  Brekhov- 
skikh’s  (1952)  criterion  for  applicability  is 

4tT/0  cos  O  »  X 

where  p  is  the  radius  of  curvature  of  the  surface  at  the  point  of  incidence 
and  ©  is  the  local  angle  of  incidence  between  the  incident  energy  vector 
and  the  unit  normal  to  the  surface  at  the  point  of  incidence . 

An  expression  for  average  backscatter  power  derived  using  the 
Kirchhoff  method  is  (Beckmann  and  Spizzichino,  1963) 

where 


£  =  /6  rr  z  ( or/  X  )Z  cos  *  O 
VK}f=  /6  7T  1  X'z  Stn  &0 

X  =  wavelength  of  incident  radiation 
cr  ~  standard  deviation  of  surface  heights 
0  -  incidence  angle  of  incident  radiation 
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L(r)  =  surface  auto -correlation  function 
T*  -  distance  on  surface 

Equation  (2.10)  is  derived  under  the  assumed  conditions  that  the 
incident  energy  is  a  plane  wave  and  the  point  of  observation  is  in  the 
Fraunhofer  zone  of  diffraction.  In  addition  the  surface  is  assumed  to  be 
perfectly  conducting,  isotropic,  and  that  no  shadowing  of  one  region  by 
another  exists.  Further,  multiple  reflection  and  specular  reflection  com¬ 
ponents  are  negligible. 

The  frequency  dependence  of  this  equation  is  determined  by  the 
argument  of  the  zero-order  Bessel  function  and  the  constant  g.  The  contri¬ 
bution  of  the  Bessel  function  to  the  integral  is  decreasing  with  increases 
in  either  angle  or  frequency.  The  factor  g  decreases  as  the  wavelength 
increases,  however  note  that  g  is  the  multiplier  in  the  exponent  which 
operates  on  the  quantity  [  1  -  C(t)J  .  Depending  on  the  value  of  a/\. 

(assuming  g»l)  the  factor  C(t)  cannot  decrease  far  from  unity  before  the 
contribution  of  the  exponent  term  becomes  negligible.  Therefore  the  region 
of  C (t)  that  contributes  is  wavelength  dependent  (Fung  and  Moore,  1966). 

This  fact  is  discussed  in  detail  in  section  5.4,  and  the  Kirchhoff  method 
is  developed  in  Appendix  A. 

2 .2  Experimental  Analyses 

Experiment  programs  designed  to  determine  the  frequency  dependence 
of  backscatter  are  few.  The  most  significant  programs  are  those  conducted 
by  Ohio  State  University  (Taylor,  1 9 5 S ;  Cosgriff,  et  al . ,  1960)  and  the  Naval 
Research  Laboratories  (Ament,  et  al.  ,  1  958;  Grant  and  Yapplee,  1957);  Sandia 
Corp.  (Edison,  et  al,  1960)  and  Johns  Hopkins  University  (Wiltse,  etal.,  1957). 
Additional  measurements  include  those  of  the  moon  (Evans  and  Pettengill,  1963; 
Evans,  1962,  1965;  Lynn,  et  al.  1964),  and  acoustic  multi-frequency  measure¬ 
ments  of  Fung  and  Leovaris  (1968) ,  and  Chapman  and  Scott  (1964) .  Studies 
presently  underway  by  the  Philco-Ford  Corp.  and  several  classified  studies 
in  progress  offer  more  promise  of  defining  the  frequency  dependence  of  back¬ 
scatter  than  any  previous  programs.  In  addition,  the  University  of  Kansas 
is  constructing  a  panchromatic  system  (Waite,  1968)  for  similar  investigations. 
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2.2.1  Ohio  State  University  Experiments 


The  OSU  Antenna  Laboratory  has  conducted  extensive  experimental 
programs  using  a  three -frequency  system  mounted  on  a  truck  platform.  The 
system  antennas  are  elevated  on  an  adjustable  boom.  The  system  operates 
at  10  GHz  (X-band),  15.5  GHz  (Ku-band) ,  and  35  GHz  (Ka-band). 

The  illuminated  area  is  generally  less  than  a  meter  in  diameter  and 
hence  the  experimental  targets  are  small  scale  terrain  features  such  as 
grass,  concrete  roads,  etc.  The  measurement  of  average  return  power  is 
obtained  by  moving  the  truck  over  the  terrain  during  the  recording  period. 

The  resulting  data  is  reasonably  accurate  and  repeatable,  and  the  terrain 
type  is  well  defined,  although  its  statistical  parameters  are  unknown.  The 
experimental  results  are  not  readily  adaptable  to  aircraft  altitude  measure¬ 
ments,  due  to  the  composite  nature  of  terrain  seen  at  these  altitudes,  but 
the  results  have  proven  very  valuable  as  a  guide  to  understanding  scattering 
from  "natural"  terrain.  Figure  2.3  shows  examples  of  OSU  data. 

2.2.2  The  Naval  Research  Laboratory  Experiments 

NRL  operates  a  four -frequency  radar  system  mounted  in  a  WV-2  air¬ 
craft.  At  each  frequency  the  radar  is  capable  of  transmitting  either  a 
horizontal  or  vertical  linearly  polarized  signal  and  can  receive  both  polar¬ 
izations  ,  hence  the  full  polarization  matrix  is  recorded.  The  system  operates 
at  428  MHz  (P-band) ,  1225  MHz  (L-band) ,  4455  MHz  (C-  band),  and  8910 
MHZ  (X-band).  The  four  antennas  have  azimuth  beamwidths  of  12°  (P-band), 

5. 5 9  (L-band),  and  5°  (C-  and  X-band). 

The  system  has  been  used  extensively  to  record  backscatter  charac¬ 
teristics  from  a  variety  of  targets,  although  the  primary  interest  has  been 
sea  clutter.  Because  of  the  narrow  elevation  beamwidths  of  the  C-  and 
X-band  antennas,  the  angular  data  are  recorded  on  multiple  passes,  con¬ 
sequently  terrain  definition  is  gross,  e.g.  "New  Jersey  Woods".  Figure 
2.4  is  an  example  of  NRL  aircraft-mounted  system  results. 

NRL  also  operated  a  three  frequency  ground  based  system  in  a  brief 
experiment  (Grant  and  Yapplee,  1957).  The  system  operated  at  wavelengths 
of  3.2  cm.  1.25  cm,  and  8.6  mm.  Measurements  were  recorded  for  various 
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Figure  2.4  Naval  Research  Laboratories  Measurements  -  New  Jersey  Woods 
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ANGLE  OF  INCIDENCE 
Figure  (a)  Terrain  covered 
with  tall  dry  weeds 


ANGLE  OF  INCIDENCE 
Figure  (b)  Nonhomogeneous 
terrain 
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ANGLE  OF  INCIDENCE 
Figure  (c)  Wet  terrain 
covered  with  tall  green 
weeds  or  flags 
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ANGLE  OF  INCIDENCE 
Figure  (d)  Tree- covered 
terrain . 


Figure  2.5  NRL  ground-based  backscatter  measurements 
(after  Grant  and  Yapplee) 


terrain  types  including  trees,  marsh,  dry  weeds,  etc.  Figure  2.5  shows 
examples  oi  measurements  obtained. 

2.2.3  Waterways  Experiment  Station  Experiments 

The  U  .S.  Army  Engineer  WES  operates  a  four-frequency  radar  system 
mounted  on  a  50-ft.  radius  arch.  The  antennas  travel  across  the  arch  on 
a  movable  platform  which  allows  viewing  of  the  target  sample  over  a  105 
degree  segment.  The  system  frequencies  are  297  MHz  (P-band) ,  5879  MHz 
(C-band) ,  9375  MHz  (X-band) ,  and  35,543  MHz  (Ka-band).  The  targets 
consist  of  specially  prepared  samples,  usually  soil,  placed  at  the  center 
point  of  the  arch  such  that  the  range  from  the  antennas  to  the  target  is 
constant  for  all  angles  of  incidence.  WES  has  facilities  for  controlled 
drying  of  the  samples  and  much  of  their  work  has  been  devoted  to  determin¬ 
ation  of  the  effect  of  moisture  content  on  radar  return.  In  addition  they 
have  dealt  extensively  with  penetration  measurements. 

Figure  2.6  shows  plots  of  WES  data  obtained  from  a  wheat  sample  at 
two  stages  of  growth. 

2.2.4  Sandia  Corporation  Experiments 

These  experiments  were  conducted  using  a  two-frequency  (415  MHz 
and  3800  MHz)  narrow-pulse  airborne  radar  system  (Edison,  et  al.,1960). 
Measurements  were  recorded  at  altitudes  of  2000  to  12000  feet  over  dif¬ 
ferent  target  areas.  The  backscatter  measurements  were  restricted  to  near¬ 
vertical  incidence  (less  than  30 c  from  the  nadir),  but  the  targets  are  reasonably 
well  documented.  The  objective  of  the  program  was  to  provide  data  for  the 
altimeter  designer.  Examples  of  the  Sandia  results  are  shown  in  Figure  2.7. 

2.2.5  Tohns  Hopkins  University  Experiments 

This  program  consisted  of  backscatter  measurements  from  the  sea 
using  multi-frequency  radar  system  mounted  on  the  bow  of  a  ship  (Wiltse, 
et  al. ,  1957) .  The  system  operated  at  9.6,  15,  24,  26,  35,  38,  and  48.7 
GHz.  Most  of  the  data  are  for  9.6,  24  and  35  GHz  and  incidence  angles  of 
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0  ,10  ,20  ,  30  0 ,  45  0  ,  60  0  ,  and  75  0  .  Figure  2 .8  shows  plots  of  the 

scattering  coefficient  as  a  function  of  frequency  for  five  incidence  angles. 

Note  that  the  scattering  coefficient  varies  as  \+3*5  for  ©  =  60°  and  \+5 

for  0  =  45  in  the  region  between  24  and  35  GHz.  This  slope  exceeds  that 

of  the  data  compiled  by  Katz  (see  Section  5.3)  which  shows  a  maximum 
+  2 

of  \  variation,  also  for  sea  return.  The  shape  of  the  frequency  dependence 
curves  are  similar  to  the  results  obtained  from  the  fundamental  surface 
measured  in  this  program  (see  Section  5.1)  which  showed  a  maximum  of 
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Figure  (a)  Wheat  Sample  -  3^  in.  height 


Figure  (b)  Wheat  sample  -  14  in.  height 


Figure  2.6  Waterways  Experiment  Station  measurements 
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CHAPTER  3 


SYSTEM  INSTRUMENTATION 


The  system  used  for  this  measurements  program  was  designed  and 
built  as  a  part  of  this  research  program.  The  need  for  flexibility  heavily 
influenced  the  selection  of  subassemblies  throughout  the  instrument,  and 
in  many  cases  the  subassembly  used  has  many  times  the  actual  capability 
required  for  this  program. 

3 . 1  System  Function 

The  system  operates  in  what  is  termed  a  "slow- sweep"  mode.  In 
this  mode  the  transmitted  frequency  is  swept  from  minimum  to  maximum  in 
a  fixed  time.  The  transmitter  is  gated  on  by  a  narrow  pulse  at  a  repetition 
rate  from  10  to  10  times  faster  than  the  sweep  rate,  hence  the  transmitter 
output  consists  of  a  pulse  train  in  which  each  pulse  contains  a  different 
carrier  frequency.  This  technique  is  illustrated  in  figure  3.1 .  The  approach 
is  actually  a  multiple- discrete  frequency  method  where  several  hundred 
discrete  frequencies  are  used  and  the  spacing  is  such  that  the  output  spec¬ 
trum  of  one  overlaps  the  spectra  ot  neighboring  pulses. 

The  envelope  of  the  average  amplitudes  of  the  returns  associated 
with  pulses  tQ  through  t2is  the  average  return  power  spectrum  for  the  parti¬ 
cular  target  being  tested.  The  handling  of  the  data  necessary  to  obtain 
this  spectrum  is  complex  and  is  described  in  detail  later. 

3.2  System  Description 

The  block  diagram  of  the  system  is  shown  in  figure  3.2.  It  consists 
of  three  main  sections:  1)  transmitting,  2)  receiving,  and  3)  data  pro¬ 
cessing  . 
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Figure  3.1  Slow- weep  mode  aampllng  technique 


Figure  3.2  Block  Diagram  of  Syatem 
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Illustration  3.1  ACOUSTIC  TEST  FACILITY 
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3.2.1  Transmitting  Section 


The  transmitting  section  consists  of  a  multi-purpose  sweep  frequency 
generator,  pulse  generator,  gate,  gated  wide-band  power  amplifier,  step- 
up  pulse  transformer,  and  transmitting  transducer. 

3.2.1 .1  Sweep  frequency  generator  —  This  generator  is  a  commercial 
Wiltron  Model  61  OB  sweep  frequency  generator  with  a  lOOKHz  to  100MHz 
Model  6104  plug-in  unit  that  serves  a  multi-purpose  role  in  the  system 
operation.  The  output  signal  can  be  swept  between  variable  f,  and  f9 
frequency  settings  in  a  variable  time.  The  output  signal  is  leveled  to 
within  +  0.25  db  across  the  range.  In  addition  to  the  output  signal,  the 
generator  provides  a  blanking  pulse,  which  is  used  to  gate  off  the  pulse 
generator  during  the  generator  retrace  interval ,  and  a  ramp  signal ,  which 
is  used  as  a  frequency  reference  in  the  processor. 

In  a  typical  operation  the  sweep  might  traverse  from  350KHz  to  1 .2 
MHz  in  5.0  seconds.  The  blanking  pulse  would  be  0.2  seconds  duration 
and  the  reference  ramp  would  rise  from  0  volts  to  10  volts  in  4.8  seconds. 

3. 2. 1.2  Pulse  Generator  —  This  generator  is  a  commercial  Hewlett  Packard 
Model  21 4A  pulse  generator  which  also  serves  a  multi-purpose  role  in  the 
system  operation.  It  provides  a  train  of  pulses  used  to  drive  both  the  gate 
and  the  gated  power  amplifier.  The  pulse  train  has  variable  pulse  width, 
amplitude,  and  period.  The  pulse  train  is  gated  off  during  the  sweep 
frequency  generator  retrace  interval  by  a  blanking  pulse  from  the  Wilton 

61 0B.  The  pulse  generator  also  provides  the  trigger  pulse  for  the  Tektronix 
scope. 

In  a  typical  operation  the  pulse  rate  is  100  pulses  per  second  with 
a  pulse  width  of  0.3  milliseconds  and  amplitude  of  -10  volts. 

3.2.1 .3  Power  Amplifier  and  gate  —  The  amplifier  is  a  Coffman  Industries 
Inc.,  Model  PA-51.  The  amplifier  supplies  55  volts  RMS  leveled  to+  1  db, 
to  a  transducer  load  over  the  frequency  range  of  lOOKHz  to  3.0MHz.  The 
gain  is  approximately  30.  The  unit  consists  of  a  wide-band,  high-power, 
push-pull  amplifier,  with  a  gain  of  10,  feeding  25  feet  of  25  ohm  coaxial 
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cable.  The  cable  is  terminated  by  a  3:1  step-up  transformer  and  200  ohm 
resistor  across  the  transducer.  An  inductor  in  series  with  the  transducer 
acts  to  form  a  low- pass  filter  with  a  cut-off  frequency  of  about  5MHz, 
depending  on  the  transducer  impedance.  This  arrangement  is  shown  in 
figure  3.3. 


wide- band 
amplifier 


25  ohm  coax 


transducer 


Figure  3.3.  Power  amplifier  termination 


The  power  amplifier  is  equipped  with  an  optional  gating  capability 
that  accepts  a  negative  5-15  volt  pulse  to  reduce  the  gain  by  biasing  off 
one  stage.  The  resulting  on-off  ratio  is  approximately  30  db.  This  ratio 
is  inadequate  to  maintain  the  receiver  input  noise  at  a  sufficiently  low 
level,  therefore  a  diode  gate  preceeds  the  power  amplifier.  This  gate  has 
an  on-off  ratio  of  40  db.  The  combined  70  db  isolation  ratio  is  satisfactory 
to  allow  adequate  signal- to- noise  ratios  for  the  measurements  recorded. 


3.2.2  Receiving  Section 

The  receiving  section  consists  of  a  receiving  transducer,  cable 
driving  amplifier,  low-noise  preamplifier,  final  amplifier,  and  signal 
attenuators . 

3. 2. 2.1  Cable  Driving  Amplifier  —  The  high  capacitive  impedance  of  the 
receiving  transducer  creates  a  major  impedance  mismatch  at  the  transducer- 
cable  connection.  This  is  partially  overcome  by  inserting  an  impedance 
transforming  amplifier.  The  amplifier  is  a  direct  coupled  two- stage  emitter 
follower  having  an  input  impedance  of  25  ohms  and  an  output  impedance  of 
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approximately  92  ohms.  The  output  is  connected  through  25  feet  of  92  ohm 
cable  to  the  low-noise  preamplifier.  The  voltage  gain  is  unity  +0.2  db 
across  the  frequency  range  0.1  MHz  to  5  MHz.  The  power  gain  is  approxi¬ 
mately  25. 

3. 2. 2. 2  Low-Noise  Preamplifier —  The  unit  is  a  Coffman  Industries,  Inc. 
Model  PA-10.  The  two- stage  preamplifier  has  a  maximum  gain  of  55  db 

+  1  db  across  the  frequency  range  0.1  MHz  to  3  MHz.  The  equivalent  input 
noise  voltage  is  less  than  30fiv.  The  input  impedance  is  matched  to  the  92 
ohm  input  cable  and  the  output  impedance  is  approximately  IK  ohm. 

3. 2. 2. 3  Final  amplifier  —  The  last  stage  of  amplification  is  obtained  in 
the  preamplifier  (Type  L)  of  a  Tektronix  545  RM  oscilloscope. 

3. 2. 2. 4  Attenuators  —  The  receiving  section  utilizes  two  variable  atten¬ 
uators  to  control  the  input  signal.  An  attenuator  prior  to  the  preamplifier 
controls  the  signal  level  to  avoid  saturation  of  the  first  stage.  Another 
0-72  db  attenuator  separates  the  two  amplifier  stages  to  protect  the 
second  stage  and  enhance  the  output  signal- to- noise  ratio.  In  typical 
operation  the  second  attenuator  is  approximately  20  db,  but  never  less  than 
lOdb. 


3.2.3  Data  Processing  Section 

The  role  of  the  data  processing  section  is  to  convert  the  amplified 
video  return  signal  to  a  form  representative  of  the  average  power  received 
as  a  function  of  frequency.  This  is  accomplished  by  sampling  the  ampli¬ 
tude  of  the  detected  signal,  converting  the  sample  from  analog  to  digital 
form,  and  storing  the  digital  information  in  a  multi-channel  memory  in  a 
manner  such  that  each  storage  channel  corresponds  to  a  different  trans¬ 
mitted  frequency.  The  section  consists  of  an  envelope  detector,  sample 
gate,  pulse  stretcher,  ramp  generator,  differential  amplifier,  gated  clock, 
comparator,  and  pulse  height  analyser. 

The  output  of  the  receiving  section  is  a  display  of  the  return  signal 
on  the  CRT  of  a  Tektronix  545B  oscilloscope.  The  input  to  the  data  pro- 
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cessing  section  is  the  vertical  signal  displayed  as  obtained  at  the  VERT 
SIG  OvJT  of  the  scope.  The  dual  time  base  of  the  545B  provides  a  variable 
width  Sumpling  gate  pulse  which  can  be  positioned  as  desired  throughout 
the  return  signal  interval.  This  gate  pulse  is  used  to  drive  a  sampling 
gate  which  selects  a  portion  of  the  return  signal  after  detection.  The 
sampled  pulse  is  stretched  and  amplitude  compared  with  a  ramp  function. 
This  operation  is  shown  in  figure  3.4. 


clock  "on"  interval 


Figure  3„4 

The  comparator  produces  an  "on"  pulse  at  the  start  of  the  stretched  pulse, 
and  an  "off"  pulse  at  the  instant  the  ramp  voltage  equals  the  sampled 
pulse  voltage.  This  interval  represents  the  "on"  time  of  a  5.0  KHz  clock. 
The  clock  output  is  a  train  of  100  n  sec,  uniform  amplitude  pulses. 

The  storage  channels  mentioned  earlier  are  actually  a  part  of  a 
256  channel  pulse  height  analyzer.  This  equipment  is  designed  to  analyze 
incoming  pulses  and  separate  them  into  channels  according  to  height.  The 
analyser  in  essence  counts  the  number  of  pulses  of  particular  heights,  i.e. 
a  single  count  is  recorded  in  one  channel  for  one  pulse  of  a  particular 
amplitude.  To  utilize  this  device  in  the  data  processing  section,  such 
that  the  channels  represent  frequency  rather  than  amplitude,  the  clock 
burst,  which  contains  pulses  proportional  in  number  to  the  signal  amplitude, 
must  be  varied  in  amplitude  across  the  sweep- frequency  range.  This  oper¬ 
ation  is  performed  by  a  comparator  which  amplifies  the  clock  pulses  by  a 
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gain  set  by  a  reference  ramp  function.  The  ramp  function  is  obtained  from 
the  sweep- frequency  generator  and  varies  proportionately  with  the  frequency 
change . 

The  record  accumulated  in  the  256  storage  channels  represents  the 
average  signal  amplitude  as  a  function  of  frequency.  The  number  of  counts 
in  each  channel  is  recorded  and  subsequently  converted  by  a  computer 
program  to  a  scattering  coefficient . 

3. 2. 3.1  Sampling  Circuit  —  The  input  to  the  processing  section  is  the 
amplified  video  return  voltage  from  the  oscilloscope  vertical  deflection 
plates.  This  signal  is  detected  in  a  linear  diode  detector  and  sampled 
by  gating  through  a  commercial  solid  state  switch,  which  is  controlled  by 
the  Tektronix  545B  +A  GATE  pulse.  The  switch  provides  an  on-off  ratio  of 
greater  than  30  db  for  the  frequency  range  50KHz  to  10MHz  with  less  than 
ldb  insertion  loss. 

3.2.3. 2  Conversion  Circuit  —  The  objective  of  this  circuit  is  to  convert 
the  amplitude  of  the  sample  pulse  into  a  time  interval  in  a  proportionate 
manner.  To  accomplish  this  the  sample  pulse  is  stretched  from  the  sample 
width  (approximately  0.02  m  sec)  to  approximately  half  the  inter- pulse 
period,  i.e.  10  m  sec  for  prf  =  50  Hz.  The  amplitude  of  the  stretched 
pulse  is  equal  to  the  amplitude  of  the  sample  pulse.  The  leading  edge  of  the 
sample  pulse  also  triggers  a  ramp  generator.  The  ramp  is  compared  with 

the  stretched  pulse  by  a  differential  amplifier  which  produces  an  output 
voltage  of  -3  volts  when  the  pulse  amplitude  exceeds  the  ramp,  and  0  volts 
when  the  ramp  voltage  exceeds  the  pulse.  The  resulting  -3  volt  pulse  has 
a  duration  proportional  to  the  amplitude  of  the  sampled  return  signal.  This 
pulse  is  fed  to  a  5  KHz  clock  "on"  gate.  The  number  of  100  n  sec  clock 
pulses  produced  is  the  measure  of  the  amplitude  of  the  return  signal. 

The  clock-pulse  train  is  one  input  of  an  operational  amplifier.  The 
reference  input  to  the  amplifier  is  the  frequency  driving  ramp  from  the 
sweep  frequency  generator.  The  ramp  duration  is  several  times  that  of 
the  pulse  train,  hence  the  output  of  the  operational  amplifier  consists  of 
a  series  of  individually  equal  amplitude  pulse  trains,  each  group  having 
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Illustration  3.2  VIEW  OF  ACOUSTIC  TEST  FACILITY 

Showing  Pulse-Height  Analyzer  (left) 


29 


an  amplitude  proportional  to  the  frequency  of  the  transmitted  signal.  The 
operational  amplifier  output  is  amplified  by  a  linear  amplifier  and  fed  to 
the  pulse  height  analyser. 

3. 2. 3. 3  Pulse  Height  Analyser  —  The  analyser  used  is  the  Radiation  Counter 

Laboratories  Mark  20  Model  2603.  It  consists  of  256  channels  of  memory 

1  fi 

having  a  per  channel  capacity  of  2  counts.  Basically  it  consists  of  an 
analog- to- digital  converter,  a  simple  digital  computer  with  magnetic  core 
memory,  a  cathode- ray- tube  data  display  system,  and  a  papers  tape  data 
printei . 

3.3  Data  Handling 

The  system  output  is  a  printed  record  of  the  number  of  counts  record¬ 
ed  in  each  of  256  storage  channels.  Each  count  represents  an  increment 
of  the  total  received  signal  voltage .  Each  channel  represents  a  segment 
of  the  total  transmitted  spectrum.  A  determination  is  made  as  to  the 
correspondence  between  channel  numbers  and  frequencies.  The  channel 
counts  are  then  smoothed  by  averaging  the  counts  in  five  channels  centered 
about  the  discrete  frequency  to  be  tabulated.  Referring  to  the  system  cali¬ 
bration  curve  of  counts  versus  voltage,  provides  the  necessary  conversion 
of  the  tabulated  results  to  voltage  and  ultimately  to  average  power  received. 

In  practice  the  record  is  preserved  as  average  number  of  counts  per 
frequency  and  later  converted  to  power.  The  reason  for  this  is  due  to  the 
calibration  procedure.  Ideally,  the  frequency  sweep  from  minimum,  f ^ ,  to 

maximum,  f9,  is  linear  and  the  traverse  from  minimum  channel,  Cf  ,  to 

1 

maximum  channel,  Cf  ,  is  also  linear  and  directly  proportional.  To  insure 

r2 

this,  the  analog- to- digital  conversion  circuit  is  set  to  provide  a  count  of 
four  for  zero  received  signal  and  a  count  of  thirty-two  for  maximum  received 
signal.  After  recording  the  return  signal  from  a  target,  a  zero  input  record 
is  made.  This  record  of  counts  is  then  subtracted  from  the  target  record, 
hence  any  non-linearity  present  during  the  test  is  eliminated.  The  count 
versus  frequency  record  is  entered  into  a  computer  program  for  conversion 
to  average  power. 
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Actually,  since  the  ultimate  measurement  desired  is  the  scattering 
coefficient,  cr°,  the  conversion  is  not  to  average  received  power,  but 
instead  to  the  ratio  of  the  average  received  power  to  average  transmitted 
power.  The  relationship  used  is: 


(3.1) 


where:  ’vPj/P^  =  power  ration  (target  signal) 

<Pr/Pt>  c  -  power  ratio  (calibration  signal) 
rc  =  distance  between  calibration  transducers 


rt  =  distance  from  transmitter  to  target 

A  =  area  illuminated  by  beam 

gt  -  normalized  gain  of  transmitting  transducer 


gf  -  normalized  gain  of  receiving  transducer 


The  calibration  information  is  the  count  record  for  the  transducer 

pair  when  they  are  signted  at  each  other  at  a  distance  r  meters.  A  more 

c 

complete  discussion  of  equation  (3.1)  is  contained  in  chapter  4. 

3.4  System  Limitations 

The  system  design  objective  of  obtaining  the  measurement  of 
broad- spectrum  backscatter  required  a  technique  of  recording  a  very 
large  quantity  of  data  quickly  and  with  a  minimum  of  manual  handling. 
This  objective  necessitated  a  relaxation  in  the  absolute  accuracy  of  the 
measurements.  In  addition  the  30:1  frequency  range  sought  was  not  ob¬ 
tained  due,  primarily  to  the  limited  capability  of  the  piezoelectric  trans¬ 
ducers  . 


3.4.1  System  Accuracy 

The  fact  that  the  scattering  coefficient  is  expressed  as  a  ratio  of 
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the  target  return  and  the  calibration  return  reduced  the  complexity  of  the 
calibration  procedure,  i.e.  whether  or  not  the  received  power  was  absolutely 
correct  was  immaterial  so  long  as  both  the  target  and  calibration  measure¬ 
ments  were  in  error  by  the  same  amount.  Even  so  the  accuracy  of  the 
scattering  coefficient  measurements  is  considered  valid  only  to  within 
3  db. 

Ihe  principal  error  is  due  to  the  method  of  data  storage.  In  practice, 
approximately  200  channels  of  the  256  Channel  Analyzer  are  employed. 

The  voltage  range  required  to  traverse  these  200  channels  is  10  volts. 
Therefore,  to  perform  the  ideal  recording  where  each  clock  count  burst  is 
recorded  in  a  single  channel,  every  pulse  of  the  possible  32  pulses  must 
be  equal  in  amplitude  to  within  +  .05  volts.  The  system  capability  is 
approximately  +  .15  volts  maximum  deviation  which  occurs  for  the  full 
count  number,  i.e.  32.  The  recording  of  the  four  count  zero  input  reference 
signal  is  almost  ideal,  but  the  subsequent  subtraction  from  the  target  sig¬ 
nal  is  unavoidably  in  error.  For  this  reason  it  was  necessary  to  smooth  the 
results  by  averaging  over  five  channels. 

A  second  cause  of  error  is  due  to  the  dynamic  range  of  the  receiving 
section  and  data  processing  section  combined.  This  dynamic  range  is 
approximately  20  db  which  is  adequate  for  the  expected  18  db  of  fading. 
However,  due  to  the  narrow  frequency  response  curve  of  the  transducers, 
much  of  this  range  was  absorbed  by  gain  changes.  By  repeating  data  runs 
several  times  and  varying  the  center  position  of  the  dynamic  ranye  with 
each  run,  it  was  possible  to  compensate  for  much  of  the  gain  change, 
however  the  extreme  ends  of  the  transducer  response  skirts  for  each  run 
were  generally  regarded  as  unusable. 

To  minimize  all  error  terms,  the  data  for  each  target  was  repeated 
until  reasonable  assurance  of  data  accuracy  was  attained. 

3.4.2  System  Frequency  Range 

The  frequency  response  of  the  transmitting  and  receiving  section 
was  discussed  earlier,  however  neither  constitutes  the  critical  control  or 
the  system  frequency  range.  This  is  controlled  by  the  transmitting  and 
receiving  transducers.  No  single  transducer  pair  is  capable  of  providing 
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sufficient  gain  to  be  usable  over  the  entire  frequency  range  desired,  i.e. 

0.1  MHz  to  3.0  MHz,  therefore  it  is  necessary  to  employ  several  pairs. 

In  most  documented  acoustic  wave  experiments,  transducer  gain 
is  measured  by  transducer  pairs  and  seldom  is  the  gain  of  an  individual 
transducer  obtained.  In  this  study  an  attempt  was  made  to  determine  the 
absolute  gain  of  a  single  transducer.  To  accomplish  this,  two  type 
measurements  were  necessary. 

3. 4. 2.1  Single  Transducer  Measurement  —  The  test  setup  for  this  measure¬ 
ment  is  shown  in  figure  3.5.  A  single  transducer  is  used  as  both  the  trans¬ 
mitting  and  receiving  antennas.  The  transducer  is  oriented  so  that  the 
acoustic  wave  propagation  is  normal  to  a  0.25  inch  thick  plate  glass  target. 

The  single  transducer  measurements  were  recorded  with  the  sweep 
frequency  oscillator  set  to  a  fixed  frequency.  The  measurement  was  the 
ratio  of  viN/vqUT-  The  transducer  gain  is  calculated  from  the  equation 


Vu. 

Vo„T 


4ir(2k) 


(3.2) 


where  2h  is  the  distance  to  the  target,  \  is  the  wave  length,  and  K  is  the 
reflection  coefficient,  which  is  unknown.  Equation  (3.2)  was  developed 
from  the  equation  for  one-way  propagation  by  assuming  the  receiving 
transducer  is  at  the  image  point  of  the  transmitting  transducer  and  that 
K  =  1 .  The  correct  value  of  K  is  determined  later. 

Plate  glass  was  used  as  the  target  to  insure  a  smooth  surface  to 
all  wave  lengths  in  the  range  of  interest.  The  calculated  reflection  co¬ 
efficient  for  glass  in  water  is  0.83. 

As  is  evident  in  the  data,  using  the  plate  glass  as  a  target  created 
some  unexpected  difficulties.  The  glass  target  was  approximately  0.25 
inches  thick,  which  is  equal  to  one  wave  length  at  about  1  MHz  (velocity 
of  propagation  of  sound  in  glass  is  6  km/sec).  Assuming  the  glass  acts 
as  a  low-loss  transmission  line,  the  terminating  impedance  on  the  backside 
of  the  plate  is  reflected  to  the  front  side  and  thus  presents  a  matched 
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Figure  3.5.  SINGLE  TRANSDUCER  TEST 
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Figure  3.6.  TRANSDUCER  PAIR  MEASUREMENT 
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Illustration  3. 


impedance  termination  to  the  incident  energy.  The  matched  condition 
should  appear  at  about  0 . 5  MHz,  1.0  MHz,  1,5  MHz,  2.0  MHz,  2.5  MHz, 
and  3.0  MHz  based  on  the  thickness  of  the  glass.  The  results  confirm  the 
condition  at  these  frequencies.  Fortunately  the  nulls  were  very  sharp  for 
frequencies  above  1  MHz.  The  two  nulls  at  1  MHz  and  0.5  MHz  were 
broad,  and  hence  the  data  are  suspect  in  this  range. 

3. 4. 2. 2  Transducer  Pair  Measurements  —  The  test  set  up  for  this  measure¬ 
ment  is  shown  in  figure  3.6.  Two  transducers  were  aligned  directly  facing 
each  other,  1 .0  meter  apart.  The  measurement  recorded  was  the  ratio  of 
voltage  at  the  receiving  transducer  to  the  voltage  at  the  transmitting 
transducer,  VQUT/VIN. 

In  this  arrangement  signals  were  recordable  for  the  entire  frequency 
range  of  0.1  MHz  to  3.0  MHz  for  nearly  all  transducers.  The  gain,  either 
receiving  or  transmitting,  is  calculated  using  the  equation  for  one-way 
propagation . 


C  -  V'N  (4uR*j 

*  V~T  GrX*  (3-3> 

where:  Gt  =  gain  of  transmitting  (or  receiving)  transducer 

Gr  =  gain  of  receiving  (or  transmitting)  transducer  measured 
by  single  transducer  method 
R  =  spacing  between  transducers 
.  \  =  wavelength 

Although  it  was  mentioned  that  the  entire  frequency  range  can  be 
covered  in  the  pair  measurements,  Gt  can  only  be  calculated  over  small 
segments  due  to  the  limited  number  of  values  of  Gr  (equation  3.2)  available 
for  calculation  from  the  single  transducer  measurements . 

The  data  plots  are  very  similar  in  shape  to  those  obtained  in  the 
single  transducer  measurement,  especially  above  1  MHz.  However, 
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there  is  a  difference  in  magnitude  of  14-15  db.  This  difference  is  due 
to  the  fact  that  the  reflection  coefficient,  K,  was  assumed  to  be  unity  in 
equation  (3 . 2).  The  actual  value  of  K  was  determined  using  both  the 
single  and  pairs  data. 


where  h  =  height  above  target  (single) 

R  =  spacing  (pairs) 

and  vqut  =  input  anc*  outPut  °*  receiving  transducer  (single) 
and  Vou^.  =  input  and  output  of  transmitting  transducer  (single) 
VjN  and  VQUT  =  input  and  output  of  pairs  measurement 

The  value  of  K  for  the  3.0  MHz  12°  and  6°  pairs  was  0.21 .  This  compares 
with  a  calculated  value  for  K  of  0.83.  The  difference  is  most  likely  due  to 
a  combination  of  effects  such  as  1)  difference  in  the  power  delivered  to 
the  transducer  mount  as  compared  to  that  actually  delivered  to  the  water, 

2)  difference  in  actual  properties  of  the  glass  target  and  ideal  properties 
of  glass;  density  equal  2400  kb/cm  and  velocity  of  propagation  equal 
6  km/sec,  and  3)  transmission  line  impedance  transformation  effects. 

Using  a  K  value  of  0.21  the  receiving  gain  was  computed  for  each 
transducer  used  in  the  experiment  both  as  individual  units  and  in  pairs. 

The  disk  or  piston  transducers  used  in  this  work  have  an  active 
surface  area  measuring  several  wave  lengths  across  and  produce  a  directive 
antenna  pattern.  The  far  field  pattern  for  a  flat  disk  transducer  operating 
in  an  infinite  baffle  is  given  by 


2  Jj(ka.,sm9) 
ka.  sm  § 


(3,5) 
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Illustration  3 


where 


Jj  =  Bessel  function  of  first  kind  and  first  order 

k  =  wave  number 
a  =  radius  of  flat  disk 
0  =  angle  from  surface  normal 

3.5  System  Imaging  Capability 

The  basic  design  of  the  system  concentrated  on  obtaining  back- 
scatter  power  spectra  from  homogeneous  rough  surfaces,  however  to  facili¬ 
tate  application  studies,  an  imaging  capability  was  also  included.  The 
system  will  record  a  B-scan  presentation  of  the  return  signal.  The  scan 
is  presented  on  a  Tektronix  561  scope  equipped  with  a  P-16  phosphor  CRT. 
This  phosphor  has  excellent  chracteri sties  for  image  recording,  including 
rapid  decay  time  and  high  dynamic  range.  The  recording  sheet  film  and  film 
developing  process  were  chosen  to  give  maximum  gray  scale  in  the  final 
image. 

The  system  arrangement  for  imaging  uses  the  same  transmitting  and 
receiving  sections,  however  the  data  processing  section  is  disconnected. 
The  VERT  SIG  OUT  of  the  Tektronix  545B  is  fed  directly  through  a  5  MHz 
band  width  amplifier  to  the  grid  of  the  Tektronix  561  CRT.  A  delayed  sweep 
time  base  unit  (3B3)  allows  positioning  of  the  return  signal  so  that  the 
intensity  modulated  raster  shows  only  the  target  return  interval.  The 
presentation  technique  employed,  as  well  as  the  type  images  obtained,  is 
discussed  in  Chapter  6. 

In  order  that  controlled  conditions  could  be  assured  during  recording 
of  images,  a  linear  motion  carriage  was  constructed.  This  device,  in 
essence,  simulates  aircraft  flight  except  that  its  motion  is  uniform  and 
extremely  stable.  The  carriage  velocity  can  be  varied  over  an  approximate 
50:1  range  with  a  minimum  of  approximately  0.5  cm/sec. 
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CHAPTER  4 


SCATTERING  EXPERIMENT 


The  experiment  program  undertakes  two  major  areas  of  investiga¬ 
tion.  This  chapter  and  Chapter  5  deal  with  the  measurement  of  broad- 
spectrum  backscatter  from  statistically  rough  surfaces.  The  objectives 
are  to  determine  the  nature  of  the  frequency  dependence  of  the  backscatter 
from  various  surfaces  and  to  compare  the  results  with  theory.  The  second 
area  of  investigation  is  handled  in  Chapter  6  which  establishes  the 
characteristics  of  broad -spectrum  illumination  applicable  to  imaging  radar 
systems . 

Two  surfaces  were  prepared  for  measurement  representing  a  gently 
undulating  "rough"  surface  and  a  gently  undulating  "rough"  surface  with 
small  scatterers  superimposed.  The  first  surface  presents  a  broad 
roughness  density  spectrum.  The  second  has  the  same  spectrum  with  a 
narrow  high  frequency  spectrum  added.  The  broad -spectrum  backscatter 
from  each  surface  was  measured  and  the  results  are  presented  in  Chapter  5. 

In  this  chapter  the  measurement  technique  is  discussed,  the  targets 
are  described,  some  characteristics  of  broad -spectrum  illumination  are 
shown,  and  the  acoustic-electromagnetic  analogy  is  explained. 

4.1  Measurement  Technique 

The  scattering  experiment  conducted  constitutes  a  case  of  re¬ 
radiation  from  a  surface  sufficiently  rough  that  the  re-radiated  energy  is 
composed  entirely  of  scatter  components .  This  is ,  the  amount  of  reflected 
energy  in  the  re-radiated  field  is  negligible.  Therefore  it  is  sufficient  to 
determine  only  the  scattering  cross  section.  The  scattering  cross  section 
is  measured  as  a  function  of  angle  and  frequency  by  measuring  the  re¬ 
radiated  power  which  is  related  to  the  scattering  cross  section  by  the 
radar  equation . 
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where:  P  (©.  X)  =  the  average  of  the  received  power  obtained  over 

an  ensemble  of  surfaces 

PT  =  transmitted  power 

GT(^p),  GR(^p)  =  the  gain  function  of  the  transmitting  and  receiving 
T  R  K  antennas  (transducers)  respectively.  Each  is 

assumed  to  be  circularly  symmetric. 

X  =  wavelength  of  radiation 

S  =  mean  surface  illuminated 


ct°  =  scattering  cross  section 


r^  =  range  from  mean 


surface  to  transmitter  and  receiver 


Backscatter  measurements  were  made  using  sufficiently  narrow 
beam  antennas  (transducers)  that  the  assumption  that  the  scattering 
cross  section  was  constant  over  the  illuminated  area  was  considered 
valid .  The  measurements  recorded  were  time  independent  over  the  effec¬ 
tive  illuminated  area  and  the  steady  state  was  assumed  to  exist.  This 
condition  is  referred  to  as  "beam  width  limitation"  and  consists  of  utilizing 
a  sufficiently  wide  transmitted  pulse  to  simultaneously  illuminate  the 
entire  effective  area  on  the  surface . 

With  the  assumption  of  constant  <j°  over  the  surface,  equation  (4.1) 


can  be  written  as: 

<P(e,x)>  = 


(4.2) 


The  calibration  procedure  detailed  in  Chapter  3  described  the  recording  of 
direct  one-way  power  with  the  transmitting  and  receiving  transducers 
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directed  toward  each  other.  This  measurement  is  convenient  for  normal¬ 
ization  of  the  average  received  power.  The  calibration  measurement  can 
be  written  as: 


p  =  Pn  GrjoKJo) 


where  rc  is  the  distance  separating  the  two  transducers.  The  normalized 
power  is  then: 


dA 


(4.4) 


where  gT  (^,)  and  gp  (^rR)  are  the  normalized  antenna  gain  functions 
having  a  maximum  of  unity.  The  integral  accounts  for  the  "aperture  effect" 
which  is  the  weighing  of  the  scattered  signals  by  the  antenna  gain  func¬ 
tions .  This  effect  has  been  detailed  by  Parkins  (1966) . 

Because  of  the  system  operation  employed  in  this  experiment  it 
is  convenient  to  express  equation  $.4) as  follows: 


cr0(Q,\)  - 


(4.5) 


The  "aperture  effect"  integral  is  calculated  separately  and  is  a  constant 

for  each  ©  and  \.  The  cross  section  is  then  determined  by  the  measured 

ratio  <P(©,  \)>  /PT  and  pTc/Pc»  The  transmitted  power  used  for  all 

measurements  and  calibrations  was  constant,  hence  P  =  P  ,  however 

1  : 

the  ratios  were  calculated  individually  and  then  multiplied  in  equation 
&.5).  The  reason  for  this  will  become  obvious  when  reviewing  the  data 


recording  procedure. 

As  described  in  Chapter  3 ,  the  backscattered  energy  is  recorded 
as  a  voltage  amplitude  on  a  frequency  ordinate,  i.e.  a  voltage  spectrum. 
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The  record  is  collected  for  the  calibration  set  up#  where  the  transmitting 
and  receiving  transducers  are  facing  each  other,  and  for  several  incidence 
angles  on  the  target.  If  the  transducer  performance  was  not  highly  fre¬ 
quency  dependent  and  their  response  were  more  repeatable,  the  target  return 
could  be  used  directly  in  equation  (4 .  ^  with  a  suitable  normalization 
factor.  However,  to  insure  valid  results  the  transducer  characteristics 
are  recorded  for  the  exact  frequency  spread  as  the  target  data  and  is  re¬ 
peated  during  each  data  run.  The  PTc/Pc  Pl<*  *or  each  data  mn  are  used 
on  a  point-by -point  frequency  basis  to  calculate  a*  in  equation  (4 .  $.  Even 
with  this  care,  the  frequency  response  curve  of  the  transducers  is  not 
sufficiently  reliable  and  repeatable  on  the  extreme  skirts  to  lend  confi¬ 
dence  to  measurements  far  removed  from  resonance.  To  illustrate  the 
characteristics  of  the  transducer,  the  pulse  height  analyzer  display  is 
shown  in  figure  4.1a  for  the  0.5  MHz  transducer  pair.  Figures  4.1b 
and  c  are  shown  to  illustrate  the  method  used  to  record  the  return  over  the 
entire  frequency  range.  The  gain  in  4.1b  is  6db  less  than  that  of  4.1a, 

4.1c  is  12db  less  than  4.1a.  Figure  4. Id  shows  the  zero  return  reference 
count.  The  spikes  at  each  end  are  the  minimum  and  maximum  frequency 
marks,  ^  and  f 2-  The  <P(©,  X)>  /PT  curve  is  calculated  for  all  three 
data  runs  and  superimposed  as  shown  in  figure  4.2  to  obtain  the  points 
used  to  calculate  the  scattering  coefficients. 

4.2  Measurement  Set-up  and  System  Parameter. 

The  measurement  of  a# ,  as  discussed  earlier,  requires  the  mea¬ 
surement  of  the  average  of  the  received  power  over  an  ensemble  of  surfaces . 
These  measurements  were  obtained  using  a  m9c.ianical  system  which 
rotated  the  statistically  rough  targets  under  the  fixed  transducers.  The 
transducers  were  mounted  on  a  boom  which  allowed  positioning  of  the 
transducers  at  fixed  angles,  relative  to  the  target  normal,  at  a  constant 

range  from  the  target. 

The  target  was  rotated  at  a  constant  rate  of  3  rpm  to  minimize  tur¬ 
bulence  of  the  water.  The  transmitter  pulse  width  of  300  microseconds 
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BACKSCATTER  MEASUREMENTS 


number  of  counts 


Figure  4.  1  Signal  return  spectra  recorded  at  three  different  receiver 
gain  settings  (a-c)  and  reference  signal  record  (d) . 
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Pr/Pt  (db) 


was  sufficiently  wide  to  give  the  desired  steady  state  condition  at  the 
extreme  angle  of  60°  incidence.  The  transmitter  PRF  was  set  a  50  Hz. 

The  sweep  rate  was  5.0  seconds,  hence  the  sample  rate  for  each 
frequency  was  0.2  Hz,  or  approximately  four  samples  per  each  rotation 
of  the  target,  which  was  sufficiently  low  to  produce  independent  data 
points.  The  data  were  averaged  for  300  seconds,  i.e.  60  sweeps. 

4 . 3  Description  of  the  Rough  Surfaces 

The  scattering  experiment  utilized  two  targets.  The  first  consis¬ 
ted  of  a  stationary  surface  of  roughness  that  was  describable  mathemati¬ 
cally  and  smooth  enough  that  locally  flat  approximation  would  be  valid. 

A  second  target  was  obtained  by  coating  the  first  target  with  uniform  dia¬ 
meter  sand  particles.  This  surface  represents  a  high  frequency,  but  dis¬ 
continuous  ,  surface  component  superimposed  on  a  low  frequency  surface 
component.  No  attempt  was  made  to  describe  this  surface  statistically . 
The  objective  was  to  introduce  a  reasonably  sharp  spike  in  the  roughness 
density  spectrum  of  a  known  surface. 

The  fundamental  surface  was  constructed  of  depolymerized  rubber 
(DPR  Compound  #242 ,  DPR  Corp.) .  This  material  was  found  to  form  a 
desirable  structural  rigidity  and  insured  minimum  propagation  of  shear 
waves.  Tests  on  this  material  indicated  its  reflection  coefficient  to  be 
constant  over  the  range  of  frequencies  used.  The  sand  subsequently  used 
on  the  surface  was  screened  to  particle  sizes  in  range  0.59 5-0. 825mm 
diameter. 

The  original  surface  was  molded  from  a  steel  surface  previously 
used  by  Parkins  (1966)  in  an  extensive  monochromatic,  bi-static  experi¬ 
ment  program.  He  determined  the  properties  of  the  surfaces  by  making 
estimates  of  the  probr.bil.ity  function  of  the  surface  heights,  the  auto¬ 
correlation  function  ol  uio  surface  heights  as  a  function  of  lateral  dis¬ 
tances,  and  the  stationarity  of  the  processes.  The  estimates  were  com¬ 
puted  from  a  series  of  sampled  height  measurements  that  were  taken  along 
profiles  of  the  target. 
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The  sampled  height  measurements  were  made  along  several  profiles 
of  the  steel  target  mold.  The  profiles  were  separated  a  distance  sufficient 
(several  correlation  distances)  to  insure  independence  of  the  data  points 
of  separate  measurement  runs .  The  probability  distribution  function  and 
auto-correlation  function  were  computed  for  each  of  these  runs  as  well  as 
the  mean#  m,  and  standard  deviation  of  heights,  a . 

The  estimate  of  the  correlation  function  was  calculated  using  the 
statistical  estimator  (the  overbar  indicates  sample  function) 


Cm 


where: 


N‘JL 


*  C  (Xax) 


th 

h^  is  the  1  height  in  the  series  of  N  points. 
X  is  the  lag . 

Ax  is  the  sampling  interval. 


h^  is  the  sample  mean. 


(4.6) 


2 

The  denominator  of  equation  (4.6)  is  the  sample  variance,  <j  ,  or  the  square 
of  the  sample  standard  deviation.  The  maximum  lag  used  in  the  calculation 
of  C(t)  was  N/10  where  N  was  made  sufficiently  large  for  the  particular 
sampling  interval  used  to  adequately  define  C(t).  The  criterion  for  the  selec¬ 
tion  of  the  sampling  interval  is  that  it  must  not  exceed  half  the  period  of 
the  highest  frequency  present  in  the  record  of  the  profile.  For  the  surface 
used  in  the  experiment,  the  value  of  Ax  =  0.02  inches  amply  satisfied 
this  criterion.  The  number  of  points  taken  in  the  measurement  was  1500 
which  corresponds  to  a  record  length  of  30  inches .  This  provided  a  maxi¬ 
mum  lag  of  one  inch  which  is  sufficient  to  adequately  define  the  falls  of 
the  sample  auto-correlation  coefficients .  The  sample  auto -correlation 
coefficients  computed  are  shown  together  with  the  sample  variances  of 
the  surface  in  figure  4.4. 
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Illustration  4.2  RUBBER  TARGET  -  FUNDAMENTAL  SURFACE 

{3  ft.  diameter) 


Illustration  4.3  HEIGHT  GUAGE  TOR  MEASURING 

SURFACE  PROFILE 
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The  behavior  of  the  sample  auto-correlation  coefficients  indicate, 
by  the  smooth  falls  and  following  rises,  that  the  surface  is  gently  undu¬ 
lating  and,  to  a  degree,  periodic.  It  was  found  that  near  the  origin  the 
sample  auto-correlation  coefficier  ^  is  closely  approximated  by  the  function 


(4.7) 


although  this  is  not  shown  in  the  figure.  However,  away  from  the  origin,  the 
curve  falls  off  more  rapidly  than  this  and  a  better  overall  approximation  is 
the  Gaussian  function 


(4.8) 


The  probability  distribution  function  of  the  sampled  heights  were 
determined  for  each  of  the  profiles  for  which  a  sample  auto-correlation 
function  was  calculated.  These  results  are  shown  in  figures  4.5  and  4.6 
together  with  plots  of  the  Gaussian  distribution  function 


2  — 

obtained  using  for  m  and  a  values  of  the  sample  mean,  m,  and  the 

_ 2 

sample  variance,  a  .  Comparison  of  the  experimental  curves  to  those 
calculated  shows  the  near-Gaussian  nature  of  the  random  height  process. 
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Run  #  1:  a  =  0.587  mm 


Figure  4.5  SAMPLE  PROBABILITY  DISTRIBUTION  FUNCTION  FOR  FUNDAMENTAL  SURFACE  RUN  #  1 

(Parkins  1965) 
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Figure  4.6  SAMPLE  PROBABILITY  DISTRIBUTION  FUNCTION  FOR  FUNDAMENTAL  SURFACE  RUN  #  2 

(Parkins  1965) 


4.4  Fading  and  Broad- Spectrum  Illumination 


The  reradiation  pattern  of  rough  surfaces  is  a  complex  arrangement 
of  maxima  and  minima  that  is  a  function  of,  among  other  things,  angle  of 
incidence  of  the  illuminating  energy.  This  complex  pattern  accounts  for 
the  phenomena  termed  fading.  For  this  discussion,  this  will  be  specified 
as  positional  fading,  that  is,  as  the  illuminating  source  is  moved  relative 
to  a  stationary  target  the  return  amplitude  will  fluxuate  pulse- to- pulse. 

The  reradiation  is  also  a  function  of  frequency,  and  definite  reorientation 
of  the  reradiation  lobes  occurs  as  the  frequency  is  changed.  This  experi¬ 
ment  attempts  to  determine  the  relationship  between  positional  fading  and 
frequency  fading . 

4.4.1  Fading  Experiment 

This  experimen  consisted  of  recording  the  return  from  the  roughened 
concrete  floor  of  the  water  tank  using  a  1 .5 MHz  transducer  pair  positioned 
at  a  45°  incidence  angle.  The  transducers  were  mounted  on  the  linear  motion 
carriage  and  pok  ted  broadside  to  the  direction  of  travel.  The  return  ampli¬ 
tudes  were  recorded  by  two  methods:  1)  the  A-scan  display  was  photographed 
on  35mm  movie  film  on  a  pulse- to- pulse  basis,  and  2 )  the  return  amplitude 
was  sampled  and  fed  to  the  pulse- height  analyzer  which  produced  the  pro¬ 
bability  density  distribution  of  the  amplitudes. 

A  position  fading  record  was  recorded  and  the  probability  density 
function  obtained.  This  distribution  is  shown  in  Figure  4.7  and  was 
obtained  by  recording  the  monochromatic  ret.’  as  for  200  seconds  as  the 
carriage  moved  the  transducers  relative  to  the  surface  at  a  velocity  sufficient 
to  assure  pulse- to- pulse  independence.  The  distribution  shown  in  4.7a 
is  clearly  Rayleigh  as  expecte  ' .  The  Rayleigh  probability  density  function 
is 


where  v  =  voltage 
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In  4.7b-e  the  system  receiver  gain  was  increased  to  observe  the  effect 
of  increased  (v  )ave  and  to  confirm  that  the  fading  was  definitely  Rayleigh. 

A  frequency  fading  record  was  recorded  and  the  probability  density 
function  obtained.  This  distribution  is  shown  in  Figure  4.8  for  four  ranges 
of  frequency  sweep.  The  carriage  is  stationary  for  this  measurement.  In 
general  no  indication  is  given  that  the  distribution  is  Rayleigh  or  Gaussian. 
The  Rayleigh  nature  of  the  distribution  shown  in  4.8d  is  misleading.  Due 
to  the  frequency  response  curve  of  the  transducers,  the  average  return 
decreases  as  the  frequency  sweep  range  is  increased,  hence,  a  nonlinear 
weighing  is  introduced  that  increases  the  predominance  of  low  amplitude 
signals.  In  figure  4.9a-c  are  shown  three  distributions  each  obtained  with 
the  same  frequency  sweep  range  but  at  three  different  set  positions  (approxi¬ 
mately  1  inch  apart)  along  the  track.  These  distributions  verify  that  the 
surface  is  not  uniformly  rough  (this  is  confirmed  by  observation  of  the 
trowel  markings  on  the  concrete  surface). 

A  position-plus -frequency  fading  record  was  obtained  by  using  a 
frequency  swept  transmitted  signal  and  recording  the  amplitudes  for  200 
seconds  as  the  carriage  moved  relative  to  the  surface.  This  record  is  shown 
in  Figure  4.10b.  and  compared  with  the  distribution  for  positional  fading 
alone  (4.10a.).  There  is  essentially  no  change  due  to  adding  frequency 
diversity.  This  result  was  confirmed  by  employing  several  carriage  veloci-. 
ties,  pulse-rate  frequencies,  and  sweep  widths .  That  is,  the  probability 
density  distribution  of  the  pulse -to-pulse  fading  of  return  from  an  area 
extensive  rough  surface  is  unchanged  by  the  fact  that  each  pulse  is  at  a 
different  frequency,  at  least  for  reasonable  frequency  ranges,  e.g.  +  20 
per  cent  from  center.  Although  it  could  not  be  confirmed  it  is  believed 
that  this  contention  will  hold  for  even  much  broader  spectrum  illumination. 

However,  even  though  the  pulse -to-pulse  fading  distribution  is 
unaffected,  the  scan-to-scan  distribution  is  greatly  altered.  That  Is  if 
the  average  return  for  each  scan  is  compared  with  the  monochromatic  pulse- 
to-pulse  return  it  would  be  seen  that  their  distributions  are  considerably 
deferent  and  that  the  scan-to-scan  fading  distribution  has  a  much  smaller 
variance.  Unfortunately,  the  system  does  not  include  the  capability  of 
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Figure  4.7  Positional  fading  probability  distributions  recorded 
at  different  receiver  gain  settings  . 


57 


pulse  counts 


a 


t  \ 


Figure  4. 


b 


voltage 


Frequency  fading  probability  distributions  for  four  ranges 
of  transmitted  spectrum  width. 
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Figure  4 
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.  9  Frequency  fading  probability  distributions  recorded  at 
different  positions. 
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Figure  4.10  Positional  fading  probability  distribution  (a)  and  frequency- 
plus -positional  fading  probability  distribution  (b). 
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automatically  producing  the  scan- to- scan  fading  distribution.  To  verify 
the  contention,  the  return  was  photographically  recorded  and  the  averages 
computed  manually.  The  sweep  time  was  set  to  be  less  than  the  interpulse 
period  of  the  monochromatic  record,  hence  the  time  between  pulses  of  the 
same  frequency  is  the  same  in  both  the  monochromatic  and  sweep  frequency 
records.  The  pulse- to-pulse  monochromatic  fading  records  had  a  mean  of 
6.83  volts  and  a  variance  of  2.967  volts.  The  scan- to- scan  swept  frequency 
(with  moving  carriage)  had  a  mean  of  6.48  volts  and  a  variance  of  only 
.545  volts. 

4 . 5  Acoustic-Electromagnetic  Analogy 

Acoustic  model  experiments  have  been  used  extensively  to  study 
a  variety  of  electromagnetic  problems .  The  validity  of  such  modelling 
experiments  depends  on  the  validity  of  the  analogy  of  the  acoustic  model 
to  the  electromagnetic  model.  The  bounds  on  analogue  will  be  discussed 
in  general  terms.  A  complete  description  of  acoustic  modeling  techniques 
used  for  radar  backscatter  is  contained  in  the  dissertation  by  Edison  (1961).  A 
detailed  examination  of  the  acoustic- electromagnetic  analogy  for  scattering 
from  rough  surfaces  meeting  the  conditions  for  application  of  the  Kirchhoff 
approximation  is  presented  in  the  dissertation  by  Parkins  (1966). 

The  analogy  between  electromagnetic  waves  in  air  and  acoustic 
waves  in  water,  is  easily  seen  by  examining  the  differential  equations  in 
each  case.  The  "telegraphers  equations"  describing  acoustic  and  electro¬ 
magnetic  plane  waves  traveling  in  the  z  direction  are: 


Acoustic  Electromagnetic 
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where: 


p  =  the  acoustic  pressure 
jj,z  =  the  particle  velocity  in  the  z  direction 
v  =  the  mass  density  of  the  acoustic  medium  per  unit  volume 
K  =  the  compressibility  of  the  acoustic  medium 
H  =  the  magnetic  field 
E  =  the  electric  field 

(x  =  the  permeability  of  the  electromagnetic  medium 
e  =  the  permittivity  of  the  electromagnetic  medium 


The  electromagnetic  equations  are  recognized  simply  as  expressions 
of  Faraday's  Law  and  Ampere's  Law  respectively.  Combining  the  two  equa¬ 
tions  in  each  case  yields  wave  equations  of  the  form 


Acoustic 


c)>.  »  L 


Electromagnetic 

cHrx  '  df* 

c)**  ’  Ci  dt* 
<9%  .  L 


where:  c  =  the  velocity  of  the  acoustic  wave  in  water 
w 

c  =  the  velocity  of  the  electromagnetic  wave  in  air 
21 

Since  identical  equations  are  obeyed  by  both  the  acoustic  and 
electromagnetic  waves,  the  solutions  are  likewise  identical  and  the  ana¬ 
logy  is  valid  for  plane  waves  in  their  respective  media. 

The  boundary  conditions  in  each  of  the  two  systems  determine  if 
the  analogy  may  be  extended  to  cover  the  reflection  and  scattering  of  the 
waves  at  an  interface.  The  boundary  conditions  are: 


Acoustic 


Electromagnetic 


■/ 


<5.  ■  4 
/4/  *  ^ 
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where  n  refers  to  the  normal  component  and  t  refers  to  the  tangential 
component;  1  and  2  refer  to  media  on  opposite  sides  of  an  interface.  The 
electromagnetic  expressions  show  the  familiar  condition  that  the  tangen¬ 
tial  components  of  E  and  H  are  continuous  across  a  boundary.  The  acoustic 
portions  express  the  less  familiar,  but  reasonable  condition,  that  in  an 
ideal  fluid,  the  pressure  on  both  sides  of  a  boundary  must  be  equal  and 
that  the  normal  component  of  velocity  must  be  continuous  across  the 
boundary. 

Only  the  scaler  form  of  the  electromagnetic  wave  equation  has 
been  used.  This  is  in  fact  the  only  form  of  the  equation  that  may  be  used 
for  the  analogy,  for  the  acoustic  equation  must  be  scaler  due  to  the  scaler 
nature  of  pressure.  The  limitations  of  the  analogy  due  to  this  fact  may  be 
seen  from  further  examination  of  the  boundary  conditions.  If  it  is  assumed 
that  the  electromagnetic  wave  is  linearly  polarized  one  may  choose  the 
analogy  such  that  the  boundary  conditions  are  the  same  in  either  system. 
Thus  the  equations  are  correct  for  horizontal  polarization  with  pressure 
analogous  with  electric  field  and  velocity  analogous  with  magnetic  field. 

By  interchanging  the  roles  of  E  and  H  in  the  equations  an  analogy  of  pres¬ 
sure  with  magnetic  field  and  velocity  with  electric  field  is  obtained.  The 
boundary  conditions  for  this  analogy  will  hold  only  for  vertical  polarization. 
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CHAPTER  5 


RESULTS  OF  SCATTERING  EXPERIMENT 


The  scattering  experiment  consisted  of  the  measurement  of  back- 
scatter  from  two  statistically  rough  surfaces.  The  fundamental  surface 
was  gently  undulating  surface  having  slopes  large  relative  to  the  incident 
wavelength.  The  other  surface  was  produced  by  adding  sand  particles  to 
the  fundamental  surface.  As  shown  by  the  results,  the  effect  of  super¬ 
imposing  small  scatterers  onto  the  fundamental  surface  was  to  convert 
the  frequency  dependent  scattering  coefficient  curve  from  a  negative 
slope  to  a  positive  slope.  An  attempt  was  made  to  explain  this  and  other 
behavioral  characteristics  and  to  relate  the  results  to  theoretical  predic¬ 
tions  of  frequency  dependence . 

The  experiment  results  presented  in  this  chapter  were  recorded  as 
described  in  Chapter  4  and,  to  the  degree  possible  by  repetitive  measure¬ 
ments,  are  considered  to  be  accurate  to  within  equipment  tolerances. 

Data  that  could  not  be  repeated  were  not  included  in  the  results.  These 
include  a  sharply  rising  behavior  of  the  scattering  coefficient  curve  in 
the  region  2.5  MHz  to  3.0  MHz  discovered  in  the  measurement  of  the 
fundamental  surface,  which  could  not  be  checked  due  to  subsequent  mal¬ 
function  of  the  transducer  pair  used  in  this  frequency  region.  This  mal¬ 
function,  and  the  unreliable  behavior  of  the  0.3  MHz  to  0 . 1  MHz  region 
of  the  low-frequency  transducer  pair,  restricted  the  range  of  the  available 
data . 

5 * 1  Measurement  Results  of  Fundamental  Surface 

The  behavior  of  the  scattering  coefficient  as  a  function  of  frequency 
for  the  fundamental,  gently  undulating  rough  surface  is  shown  in  figures 
5.1  through  5.5  at  incidence  angles  of  0°  ,  15°  ,  30c  ,  45°  and  60°  . 

Figure  5.6  shows  all  the  curves  plotted  together. 

The  curves  show  two  behavioral  regions;  one  below  and  one  above 
about  1.25  MHz  (\  =  1.2mm).  In  the  region  below  1.25  MHz  the  0°  and 


15°  data  have  a  negative  slope  which  increases  above  1.25  MHz.  For 
30°  ,  45°  ,  and  60°  the  lower  frequency  region,  at  least  down  to  0.5  MHz 
(X  =  3.0mm) ,  the  slope  is  approximately  zero.  Below  0.5  MHz  a  positive 
slope  is  suggested,  especially  at  45°  and  60°  ,  however  the  confidence 
level  of  these  data  is  not  as  high  as  that  above  0.5  MHz. 

The  region  above  1.25  MHz  has  a  pronounced  negative  slope  for 
all  angles.  This  slope  is  confirmed  to  2.5  MHz  (\  =  0.6mm)  and  is 
believed  to  continue  to  higher  frequencies,  however  some  measurements 
indicated  a  rapid  change  to  a  large  positive  slope  in  the  region  above 
2.5  MHz.  These  measurements  could  not  be  confirmed  and  no  physical 
explanation  was  found  for  such  a  behavior,  therefore  the  data  were  dis¬ 
carded  . 

The  data  were  also  plotted  as  a  function  of  angle  for  several  dis¬ 
crete  frequencies  in  figure  5.7.  The  plots  indicate  a  greater  negative 
slope  near  vertical  for  the  low  frequencies  with  a  gradual  flattening  out 
to  about  1.2  MHz  (X=  1.2mm).  Thereafter  the  plots  hold  approximately 
the  same  shape  but  gradually  drop  in  overall  amplitude. 

5.2  Measurement  Results  of  Sanded  Surface 

The  frequency  dependence  of  the  scattering  coefficient  measured 
for  the  fundamental  surface  densely  coated  with  0 .595-0 . 825  mm  sand 
particles  are  shown  in  figures  5.8  through  5.12  at  incidence  angles  of 
0°  ,  15°  ,  30°  ,  45°  ,  and  60°  .  Figure  5. 13  shows  all  curves  plotted  on 
the  same  graph. 

As  in  the  case  of  the  fundamental  surface,  these  data  exhibit  dif¬ 
ferent  behavior  in  two  frequency  regions  .  In  the  region  below  1.25  MHz 
the  curves  have  a  positive  slope  of  approximately  the  same  magnitude 
for  all  angles.  Above  1.25  MHz  the  data  tend  to  become  frequency  inde¬ 
pendent.  There  was  some  support  for  defining  a  third  region  above  about 
2.0  MHz  (X  =  0.75  mm)  where  the  data  exhibit  a  positive  slope  such  as 
was  found  for  the  fundamental  surface.  Unfortunately  this  could  not  be 
confirmed  due  to  the  malfunction  of  the  transducer  pair  used  in  this  region. 
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FREQUENCY  DEPENDENCE  O l  SCATTERING  COEFFICIENT 
0  Degree  Incidence  Angle 
Fundamental  Surface 
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FREQUENCY  DEPENDENCE  OF  SCATTERING  COEFFICIENT 
30  Degree  Incidence  Angie 
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FREQUENCY  DEPENDENCE  OF  SCATTERING  COEFFICIENT 
60  Denree  Incidence  Angle 
Fundamental  Surface 
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Figure  5.7  Variation  of  scattering  coefficient  with  angle  and 
frequency  -  fundamental  surface. 
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FREQUENCY  DEPENDENCE  OF  SCATTERING  COEFFICIENT 
0  Degree  Incidence  Angle 
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FREQUENCY  DEPENDENCE  OF  SCATTERING  COEFFICIENT 
1 S  Degree  Incidence  Angle 
Fundamental  Surface 
(Covered  with  Sand  Particles) 
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FREQUENCY  DEPENDENCE  OF  SCATTERING  COEFFICIENT 
30  Degree  Incidence  Angle 
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FREQUENCY  DEPENDENCE  OF  SCATTERING  COEFFICIENT 


77 


Frequency  (MHz) 


FREQUENCY  DEPENDENCE  CF  SCATTERING  COEFFICIENT 
Fundamental  Surface 
(Covered  with  Sand  Particles) 


Frequency  (MHz) 


Scattering  Coefficient  (db) 


These  data  were  also  plotted  as  a  function  of  angle  for  several 
discrete  frequencies  in  figure  5.14.  The  plots  indicate  a  greater  nega¬ 
tive  slope  near  vertical  for  the  low  frequencies  with  a  gradual  flattening 
as  the  frequency  increases  out  to  about  1.25  MHz  (\  =  1.2  mm)  just  as 
in  the  case  of  the  fundamental  surface.  The  measurements  at  45°  inci¬ 
dence  confuse  the  interpretation  of  the  data  and  an  immediate  explanation 
is  not  available  for  this  phenomenon.  The  readings  tend  to  conform  to 
expectations  as  the  frequency  increases.  Above  1 .25  MHz  the  curves 
remain  almost  unchanged  within  the  equipment  accuracy  tolerances . 

5.3  Comparison  of  Results  with  Experiments 

Following  the  analysis  procedure  used  extensively  by  Katz  (1966) 
the  scath~  Ing  coefficient  frequency  dependence  curves  were  fit  to  a  func¬ 
tion:  a0  =  a\a.  The  fit  was  done  by  frequency  regions  and  a  determination 
of  the  best  a  was  made  for  all  data.  In  figures  5.15  and  5.16  the  variation 
of  the  coefficient  a  is  plotted  as  a  function  of  incidence  angle  for  the  funda¬ 
mental  surface.  In  figures  5.17  and  5.18  the  coefficients  are  plotted 
for  the  sanded  surfaces. 

In  the  next  section  an  attempt  will  be  made  to  explain  the  behavior 
of  the  curves  using  the  guidelines  offered  by  the  Kirchhoff  theoretical 
method  of  expressing  the  scattering  processes,  however  as  reference  to 
other  experimental  results,  the  values  of  the  coefficient  for  angular 
backscatter  measurements  of  sea  clutter,  cities,  snow  covered  terrain, 
grass  surfaces,  the  moon,  and  asphalt  and  concrete  roads  are  shown  in 
figure  5.19.  These  data  were  compiled  by  Katz  (1966)  using  data  from 
the  Ohio  State  University  and  Naval  Research  Laboratory  measurements 
programs.  These  data  were  obtained  by  polychromatic  systems  (refer 
to  Chapter  2) ,  hence  the  slope  factors  were  obtained  using  widely  spaced 
points.  Such  subtleties  as  were  found  in  the  measurements  of  the 
surfaces  used  in  this  program  could  not  have  been  detected. 
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Figure  5.16  Variation  of  slope  for  fundamental 
surface  data  above  1.25  MHz 
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Figure  5.17  Variation  of  slope  for  sanded  surface  data 
below  1.25  MHz 


Figure  5.18  Variation  of  slope  for  sanded  surface  data 
above  1.2  5  MHz 
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(a)  Wavelength  dependence  of  grass, 
concrete,  and  asphalt  surfaces. 
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(b)  Wavelength  dependence  of  sea 
clutter.  This  figure  covers  the 
microwave  region  from  0.86  to 
71  cm  and  applies  to  rough  sea 
conditions . 
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(c)  Wavelength  dependence  of  the 
moon's  surface. 


Figure  5.19  Variation  of  slopes  for  natural  surfaces 
(after  Katz  1966) 
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5.4  Comparison  of  Results  with  Theory 


Several  authors  have  observed  that  the  scattering  behavior  of 
rough  surfaces  varies  with  wavelength,  but  a  general  expression  for  this 
variation  is  not  available.  Various  theories  have  indicated  a  wavelength 

i  2  m m  c 

variation  of  \  to  X  (Fung  and  Moore,  1964;  Beckman,  1965).  Terrain 

+  5  -14 

measurements  have  shown  variations  of  X  to  X  (Katz,  1966;  Janza, 

» 

1963;  Wiltse,  et  al.  1957).  The  results  presented  for  the  fundamental 

surface  (Figure  5.15  and  5.16)  show  variations  of  X  to  X  ,  and  X 
_  2 

to  X  for  the  sanded  surface  (Figure  5.17  and  5.18). 

Even  though  a  general  expression  is  not  available,  it  is  normally 
accepted  that  surfaces  tend  to  be  "rougher'*  as  wavelength  decreases, 
and  the  angular  variations  of  the  scattering  coefficient  tend  to  flatten  so 
that,  for  an  extremely  "rough"  surface,  cj°  is  independent  of  angle  and 
frequency . 

The  characteristics  of  the  fundamental  surface  described  in  Section 
4.3  are  such  as  to  assume  the  Kirchhoff  method  to  be  an  applicable  means 
of  describing  the  backscatter.  In  Section  2.1  an  expression  for  the  average 
backscatter  power  based  on  the  Kirchhoff  method  was  presented  as  equation 
(2.10). 


(2.10) 


In  Appendix  A  the  Kirchhoff  method  is  outlined  and  reduced  to 
expressions  for  average  backscatter  power  for  "slightly  rough"  (equation  A. 
39)  and  "very  rough"  (equation  A. 40)  normally  distributed  surfaces  having 
Gaussian  autocorrelation.  For  the  case  where  g  is  large  relative  to  unity, 
the  surface  is  termed  "very  rough"  and  the  expression  is  (omitting  a  constant) 


(5.1) 

J 
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where  T  =  correlation  distance 
2 

a  =  variance  of  surface  heights 
©  =  incidence  angle 

Equation  (5.1)  is  independent  of  wavelength  and,  even  though  the  measured 
autocorrelation  of  the  fundamental  surface  appeared  Gaussian,  the  expression 
cannot  be  used  to  describe  the  frequency  dependence  measurements. 

In  Section  A. 5  the  Kirchhoff  method  is  applied  to  obtain  expressions 
for  average  backscatter  power  for  "slightly  rough"  (equation  A„47)  and 
"very  rough"  (equation  A .  48)  normally  distributed  surfaces  having  exponential 
autocorrelation.  For  g  large  relative  to  unity  ("very  rough"  surface)  the 
expression  is  (omitting  a  constant): 


where 


From  profile  measurements  of  the  fundamental  surface  the  standard  deviation 
was  0.71  mm  and  the  correlation  distance  was  17.8  mm  (run  #  2).  Using 
these  values,  equation  (5.2)  plots  as  shown  in  Figure  5.20.  These  plots  do 
not  fit  the  measured  results,  especially  for  incidence  angles  greater  than 
30° .  It  was  noted  in  Section  5.1  that  the  data  are  seemingly  separated  into 
different  behavior  regions  below  and  above  about  1 .25  MHz.  Since  the 
maximum  ot  predicted  by  any  of  the  above  expressions  is  +  2,  which  is  the 
slope  of  the  zero  degree  incidence  curve,  the  region  above  1 .25  MHz  cannot 
be  fit  (refer  Figure  5.16).  However  using  the  run  #  1  profile  measurements  of 
cr  =  .678  mm  and  T  =  13.97  mm  the  theoretical  curves  for  0°  and  15°  incidence 
angles  are  fit  reasonably  well  in  the  region  below  1.25  MHz,  but  the  fit  to 
30°  is  poor  (Figure  5.21).  These  values  represent  a  (T /cr  )  value  of  29.66. 
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It  was  noted  that  a  (T /a  )  value  of  20. 97  provided  a  reasonable  fit  to  the 
15°  and  30°  measurements  in  the  region  below  1.25  MHz,  but  the  fit  to 
0°  is  poor  (Figure  5.22).  Attempts  to  fit  the  data  above  30°  incidence  do 
not  appear  practical  using  equation  (5.2),  but  the  literature  shows  that 
there  is  no  general  agreement  as  to  the  validity  of  the  theory  for  large 
angles  of  incidence. 

The  fact  that  different  angular  regions  were  fit  by  use  of  different 
standard  deviations  and  correlation  distances  was  also  noted  by  Fung  and 
Moore  (1964)  in  attempting  to  fit  lunar  backscatter  data  using  equation 
(5.2).  The  fact  that  these  parameters  appear  to  be  angle  dependent  suggests 
that  they  may  also  be  wavelength  dependent.  This  has  also  been  suggested 
by  several  authors  (Fung  1965)  who  note  an  apparent  "size  filtering"  effect 
due  to  the  incident  wavelenth. 

To  examine  this  further,  it  will  be  assumed  that  the  ratio  (T/cr  )  is 
an  effective  parameter  which  is  a  function  of  angle  and  wavelength.  Further, 
it  is  arbitrarily  assumed  that  the  dependence  is  of  the  following  form: 


where:  R 

e 


2 

R  =  actual  parameter,  (T/cr  ) 

f  (©)  =  undefined  function  of  incidence  angle 

X  =  incident  wavelength 

\n  =  wavelength  for  which  R  =  Rf(0) 
u  s 


(5.3) 


If  p  is  positive,  as  intuitively  it  should  be,  then  equation  (5.2) 
predicts  a  slope  coefficient,  a,  greater  than  +2,  and  hence  the  region  above 
1.25  MHz  could  be  explained.  In  the  region  above  1 .25  MHz,  the  0°  data 
vary  as  X+2‘®.  Using  (5.3)  in  (5.2)  leads  to  an  expression  which  will  fit 
the  measurement  results  if  p  =  0.4.  Figure  5.23  shows  the  comparison  of 
the  measured  and  predicted  results  using  p  =  0.4,  Xq  =  1.25  mm,  R  =  3C..317, 
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Figure  5 . 20 

PREDICTED  FREQUENCY  DEPENDENCE 


Frequency  (MHz) 


COMPARISON  OF  RESULTS  WITH  THEORY 
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Frequency  (MHz) 


and  f(0)=  1  for  the  region  above  1 .25  MHz  (X  ■=  1.25  mm).  This  approach 
improved  the  fit  to  the  0°  data,  however  the  15’  and  30°  data  fit  remains 
poor . 

5.4.1  Summary  of  Comparison  with  Theory 

Although  some  degree  of  matching  of  the  measured  and  theoretical 
results  was  obtained,  considerable  doubt  exists  as  to  the  ability  of  the 
Kirchhoff  method  to  predict  the  frequency  dependence  of  backscatter,  at 
least  in  cases  where  the  actual  surface  autocorrelation  is  approximated  by 
an  autocorrelation  function.  A  point  of  major  concern  is  that  the  best  fit 
occurs  using  the  expressions  derived  for  an  e  xponential  autocorrelation 
function,  which  clearly  cannot  fit  the  measu/ad  surface  autocorrelation. 

This  contradiction  is  somewhat  overcome  by  the  approach  used  by  Fung  and 
Leovaris  (1968),  which  in  essence  maintains  the  form  of  equation  (5.2), 
but  allows  incorporation  of  autocorrelation  functions  other  than  the  exponential. 
This  approach  is  discussed  in  Appendix  B.  Even  with  this  approach  the 
maximum  wavelength  variation  is  X  (at  0  =  0°),  and  the  coefficient  must 
decrease  as  0  increases.  This  characteristic  is  contradicted  by  the  measure¬ 
ments  of  the  fundamental  surface  (see  Figure  5.6)  and  the  results  of  Wiltse, 
et  al.  (see  Section  2.2.5).  Likewise  the  theory  cannot  predict  results  such 
as  those  compiled  by  Katz  (see  Section  5.3). 

The  effective  parameter  approach  is  supported  by  the  physical 
reasoning  that  different  subranges  of  structure  sizes  dominate  the  return  at 
different  incidence  angles  and/or  wavelengths.  These  subranges  have 
autocorrelations  and  variances  which  are  not  the  same  as  the  autocorrelation 
and  variance  of  all  subranges  combined,  i.e.  the  entire  surface.  Several 
attempts  have  been  made  to  incorporate  this  character  into  the  scattering 
expressions  to  describe  the  angle  dependence  (Beckman  i965;  Fung  and 
Moore  1964,  1966;  Davies  1954).  These  attempts  have  been  only  moderately 
successful  and  the  use  of  equation  (5.3)  is  likewise  only  moderately 
successful. 

In  summary,  the  measurements  obtained  in  this  program,  as  well  as 
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Figure  5.23 

COMPARISON  OF  RESULTS  WITH  THEORY 
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Frequency  (MHz) 


those  of  other  programs,  suggest,  in  agreement  with  Davies  (1954),  that 
use  of  the  autocorrelation  function  approximation  of  the  surface  auto¬ 
correlation  does  not  give  sufficient  statistical  information  to  describe  the 
scattering  from  the  surface. 
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CHAPTER  6 


IMAGE  EXPERIMENT 


In  recent  years  the  use  of  imaging  radar  systems  for  both  military 
reconnaissance  and  geoscience  investigations  has  expanded  rapidly. 
Associated  with  this  expanding  demand  has  been  the  growing  realization 
that  the  monochromatic  nature  of  such  systems  handicaps  their  utility, 
especially  in  military  reconnaissance  applications.  This  chapter  presents 
an  experiment  with  a  panchromatic  imaging  system.  In  this  experiment 
monochromatic  images  of  complex  targets  are  compared  with  panchromatic 
images  to  show  the  effect  of  frequency  averaging.  The  substantial  improve¬ 
ment  in  target  definition  obtained  is  impressive,  especially  in  view  of  the 
fact  that  less  than  a  +  10%  bandwidth  is  employed,  ar.d  that  little  increase 
in  system  power  or  complexity  is  required  to  gain  this  improvement. 

The  theme  of  the  chapter  is  illustrated  in  Figure  6.1  which  shows 
an  acoustic  image  of  two  spheres.  The  images  actually  show  an  expanded 
view  of  a  small  area  of  the  sphere's  surface  near  the  apex.  The  return 
signal  strength  from  regions  away  from  the  apex  is  below  the  noise  level. 

In  6.1  (a)  the  spheres  are  imaged  monochromatically  and  lobing  is  evident. 
In  6.1  (b)  the  frequency  was  swept +10%  from  the  center  frequency.  The 
frequency  averaging  of  the  lobes  is  obvious.  The  monochromatic  lobes 
are  investigated  in  detail  in  section  6.3  and  the  effect  of  frequency  aver¬ 
aging  is  illustrated. 

6.1  Measurement  Technique 

Acoustic  images  were  obtained  in  a  manner  completely  analogous 
to  an  airborne  imaging  radar  system.  The  antennas  (transmitting  and 
receiving  transducers)  were  moved  past  the  targets  at  a  fixed  distance 
and  uniform  velocity.  The  transducers  were  mounted  on  a  linear  motion 
carriage  that  was  capable  of  a  wide  range  of  movement  velocities  and 
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(a)  Monochromatic  image  of  two  spheres 
(f=  1 .5  MHz) 


(b)  Swept  frequency  image  of  two  spheres 
(f=  1.5  MHz  +  10%) 


1  Effect  of  frequency  average  on  image  of  sphere  targets 
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transducer  position  settings.  The  position  geometry  is  shown  in  Figure 
6.2  and  a  photograph  of  the  arrangement  is  also  shown. 

The  system  operation  associated  with  obtaining  acoustic  images 
was  described  in  Section  3.5.  In  review,  the  signal  return  is  envelope 
detected  and  applied  to  the  grid  of  an  imaging  CRT.  Each  return,  i.e. 
the  return  associated  with  each  transmitted  pulse,  intensity  modulates 
one  horizontal  line  trace  on  the  CRT.  Sequential  returns  form  sequential 
intensity-modulated  line  traces  each  vertically  displaced  from  each  other. 

The  end  result  is  an -intensity  modulated  raster  display.  A  photographic 
image  is  obtained  by  exposing  the  film  during  the  entire  vertical  traverse 
time  interval.  The  procedure  is  illustrated  in  Figure  6.3. 

The  imaging  experiment  was  designed  around  the  limitation  presented 
by  the  transducers.  That  is,  the  available  piston  transducers  provide  only 
conical  beams  and  have  relatively  high  sidelobes.  The  use  of  spherical 
targets  insured  that  the  sidelobes  of  the  antenna  did  not  influence  the  return 
over  the  area  of  investigation  of  the  reradiation  pattern.  Ideally  the  exper¬ 
iment  would  concentrate  on  imaging  area-extensive  targets  such  as  those 
of  concern  to  imaging  radar  systems,  however  the  lack  of  a  fan-beam  antenna 
discouraged  this  approach.  Fortunately,  as  will  be  shown  in  Section  6.2, 
the  spherical  targets  provided  the  information  needed  to  relate  to  the  area 
extensive  imaging  problem,  hence  the  limitation  to  a  conical  beam  did  not 
restrict  the  value  of  the  experiment. 

The  parameters  of  the  experiment  were  set  with  the  following  con¬ 
sideration.  The  objective  of  the  experiment  was  to  determine  the  effect 
of  introducing  a  broad -spectrum  signal  in  place  of  a  monochromatic  signal. 
The  targets  examined  were  complex  and  consequently  provide  a  reradiation 
pattern  consisting  of  a  multiple  of  fine  lobes.  To  determine  the  effect  of 
frequency  averaging  required  that  the  effect  of  sample  averaging  be  excluded. 
That  is,  th~  image  would  be  altered  if  the  number  of  "looks"  were  changed 
by  the  fact  that  more  or  less  lobes  would  be  seen.  To  eliminate  this 
variable,  the  sampling  rate  (PRF)  was  set  sufficiently  high  that  all  lobes 
would  be  recorded  and  the  rate  was  not  altered  when  the  swept  frequency 
was  introduced.  The  carriage  velocity  was  set  at  6.5  mm/sec  and  the  PRF 
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Figure  6.2  Positioning  for  acoustic  Imaging. 
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Figure  6.3  Image  recording  of  target  return. 
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Figure  6.4  8-element  array. 
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illustration  6.2  LINEAR  MOTION  CARRIAGE 


was  set  at  approximately  130  Hz.  This  provided  one  sample  for  each 
.05  mm  of  carriage  travel. 

In  Chapter  3  the  system  operation  was  described  and  the  "slow- 
sweep"  mode  was  outlined.  This  mode  of  operation  consists  of  sampling 
t  ->  swept  frequency  at  the  pulse-rate  frequency.  Ideally  broad-spectrum 
imaging  would  be  done  using  a  "fast-sweep"  mode,  that  is  the  entire 
spectrum  would  be  transmitted  during  each  pulse.  This  could  conceivably 
be  simulated  by  increasing  the  PRF  proportional  to  the  sweep  time.  That 
is,  if  the  PRF  for  monochromatic  imaging  is  100  Hz,  set  the  sweep  time 
at  (1/100)  sec  and  increase  the  PRF  to  400  Hz  to  provide  4  samples  per 
sweep.  Unfortunately  this  reintroduces  the  effect  of  increasing  the  sampling 
rate  and  hence  for  this  experiment  the  PRF  for  both  monochromatic  and 
swept  frequency  signals  was  the  same,  130  Hz.  The  sweep  time  was 
adjusted  to  give  4  samples  per  sweep  (about  23.5  m  sec).  The  samples 
were  not  synchronized  with  the  sweep,  and  hence  the  frequency  sampling 
varies  sweep-to-sweep.  In  this  manner  the  full  sweep  is  effective  even 
though  the  sampling  rate  is  very  low. 

The  width  of  the  pulse  was  set  at  .6  -  .8  msec  to  provide  a  return 
segment  which  did  not  include  the  effect  of  the  high  frequency  components 
in  the  leading  and  trailing  edge  of  the  transmitted  pulse.  It  was  determined 
that  the  spherical  targets  caused  the  return  signal  to  be  restricted  to  an 
area  near  the  apex,  hence  the  leading  edge  effect  decayed  rapidly.  There¬ 
fore  by  extending  the  pulse  width  beyond  this  decay  interval,  the  leading 
edge  effect  was  excluded.  This  effect  is  evident  in  the  images  shown  in 
Figure  6.1. 

The  useful  bandwidth  of  the  swept  signal  was  restricted  by  the 
frequency  response  of  the  transducers.  The  response  curves  fall  so  rapidly 
that  variations  greater  than  about  +  10%  from  resonance  were  impractical. 

6.2  Targets 

Three  primary  targets  were  used  in  this  experiment:  (1)  a  4.5  cm 
diameter  rubber  sphere;  (2)  a  7.8  cm  diameter  styrofoam  sphere;  and  (3) 
an  8-element  planar  array  consisting  of  2.5  cm  diameter  wooden  spheres 
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Illustration  6.3  IMAGING  EXPERIMENT  OBJECT  TARGE iS 


seated  on  3.2  cm  centers  in  two  rows  of  4-elements  each.  The  targets 
are  shown  in  the  photograph  as  they  were  positioned  in  the  tank.  The 
rubber  sphere  is  on  the  left. 

6.2.1  Rubber  Sphere 

This  target  is  a  relatively  smooth,  soft  rubber  sphere  selected  to 
provide  a  broad,  hopefully  frequency-insensitive,  lobe  in  the  direction  of 
backscatter.  Although  the  exact  composition  of  the  sphere  is  unknown, 
it  is  reasonable  to  assume  negligible  shear  waves  and  a  reasonably  high 
reflection  coefficient.  The  diameter  of  the  sphere  is  very  large  relative 
to  the  incident  wavelength,  but  is  smaller  than  the  illuminated  area  of  the 
beam.  Under  these  conditions  the  scattering  cross  section  should  be  almost 
independent  of  wavelength  (Stratton,  1941)  and  the  sphere  would  provide 
the  necessary  calibration  reference  for  other  measurements. 

The  measurements  presented  in  Section  6.3  verify  the  presence  of 
a  strong,  almost  frequency  insensitive,  reflection  in  the  backscatter  direc¬ 
tion,  however  the  fact  that  the  sphere  is  not  perfectly  smooth  to  the  incident 
wavelength  is  evident. 

6.2.2  Styrofoam  Sphere 

This  target  was  selected  to  illustrate  an  effect  very  common  to 
imaging  systems  but  one  that  is  difficult  to  simulate  acoustically  in  a 
controlled  manner;  that  is,  an  area  extensive  surface  consisting  of  various 
size  randomly  oriented  scatterers.  The  surface  of  the  sphere  is  rough  in 
the  somewhat  facet-like  manner  typical  of  styrofoam.  Consequently  some 
scatterers  are  sufficiently  close  together  to  form  highly  directional  reradi¬ 
ation  patterns  similar  to  closely  spaced  arrays.  In  addition  the  multilobe 
pattern  of  widely  spaced  elements  exists,  and  scatterer  sizes  both  large 
and  small  relative  to  the  incident  wavelength  are  present.  The  spherical 
configuration  reduces  the  illumination  area  edge  effects  and  helps  accent 
the  behavior  of  the  surface  reradiation  lobes  by  restricting  the  contributing 
area  to  a  much  smaller  size  than  is  otherwise  possible  using  the  available 
acoustic  transducers . 
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6.2.3  Array 

The  array  target  consists  of  eight  spheres  arranged  as  shown  in 
Figure  6.4.  The  purpose  of  this  target  is  to  provide  a  large  number  of 
narrow  lobes  of  different  amplitudes.  The  array  is  an  extension  of  a 
4-element  array  (one  row  of  spheres)  which  provided  the  desired  fire 
lobe  structure,  but  the  lobes  were  of  almost  equal  emplitude.  The  effect 
of  variations  in  incident  wavelength  is  better  illustrated  by  use  of  the 
more  directional  pattern  of  the  8-elament  array. 

The  reradiation  pattern  of  an  array  of  point  scatterers  spaced  several 
wavelengths  apart  can  be  calculated,  and  the  results  show  a  number  of  lobes 
proportional  to  the  number  of  elements .  The  width  of  each  lobe  is  also  pro¬ 
portional  to  the  number  of  elements  (Silver  1964).  The  spherical  elements 
of  the  target  array  cannot  be  interpreted  as  point  scatterers;  however  the 
area  on  the  sphere  contributing  to  the  reradiation  is  sufficiently  small  that 
a  very  similar  pattern  to  that  calculated  is  created. 

6.3  Results 

The  results  clearly  indicate  that  the  addition  of  frequency  averaging 
of  the  reradiation  pattern  of  complex  targets  substantially  improves  the  de¬ 
finition  of  the  target.  The  reader  should  bear  in  mind  that  the  results  pre¬ 
sented  represent  only  the  addition  of  a  swept  frequency  during  transmission, 
done  in  such  a  manner  that  each  pulse  modulated  a  different  carrier  fre¬ 
quency  on  a  repetitive  sweep  basis.  Several  other  techniques  are  possible, 
some  perhaps  more  optimum  than  the  method  used  hers  (  Gustafson  and  As, 
1964;  Kosowsky,  et  al. ,  1963;  Ray,  1966).  However,  the  point  is  well 
made  that  the  addition  of  frequency  averaging  warrants  consideration  in 
the  design  of  imaging  radar  systems. 

6.3.1  Rubber  Sphere 

Figure  6.5  shows  the  acoustic  image  of  the  rubber  sphere  for  four 
separate  frequencies  and  a  swept  frequency.  The  images  were  recorded 
from  a  distance  of  51  cm. ,  and  consequently  the  illumination  region  of  the 
transmitted  beam  was  6  cm  in  diameter.  The  target  diameter  was  4 .5  cm. 


The  images  of  the  styrofoam  sphere  and  the  array  show  the  exis¬ 
tence  of,  what  may  properly  be  termed,  multiple  lobes.  The  return  from 
the  rubber  sphere  consists  of  multiple  "specular  points".  That  is,  the 
sphere  is  obviously  not  perfectly  smooth  to  the  incident  wavelength, 
but  the  roughness  does  not  produce  an  array  effect  such  as  the  styrofoam 
sphere.  The  significance  of  this  distinction  is  that  for  an  array  the  lobes 
"move"  with  frequency  and  hence  an  averaging  is  obtained.  This  effect 
does  not  occur  with  the  rubber  sphere.  The  distinction  evident  between 
images  at  different  single  frequencies  is  that  each  "specular  point"  has 
a  different  amplitude  from  frequency-to- frequency .  The  swept  frequency 
image  clearly  shows  the  target  to  have  two  major  "specular  points" ,  where 
one  is  6  db  greater  than  the  other.  Other  such  points  are  indicated  to 
exist  at  lower  levels.  No  single  frequency  image  clearly  defines  both 
points,  and  limited  to  any  one  single  frequency  image  the  definition  of 
the  target  is  misleading.  In  comparing  figure  6.5  to  figure  6.1,  both  of 
which  show  the  rubber  sphere,  the  reader  shouL  note  that  figure  6.1  was 
obtained  by  imaging  a  region  of  the  sphere  slightly  off-center  and  hence 
the  dynamic  range  of  the  "specular  point"  return  was  less  than  the  returns 
shown  in  figure  6.5,  which  are  from  the  apex. 

6.3.2  Styrofoam  Sphere 

Figure  6.6  shows  the  acoustic  image  of  the  styrofoam  sphere  for 
several  single  frequencies.  In  Figure  6.7  one  of  the  monochromatic  images 
is  compared  to  swept  frequency  images  of  +  3%  and  +6%  from  center  fre¬ 
quency.  The  images  were  recorded  at  a  distance  of  approximately  1  m 
and  a  beam  illumination  region  of  11. 6  cm  diameter.  The  reradiation  lobe 
pattern  is  more  complex  than  the  simple  rubber  sphere  and  the  effect  of 
varying  frequency  is  marked.  Of  special  interest  is  the  phase  cancella¬ 
tion  of  the  center  main  lobe  (well  defined  at  1520  KHz  and  162  6  KHz)  at 
1427  KHz,  1744  KHz,  and  1816  KHz.  The  effect  occurred  at  numerous 
intermediate  frequencies  as  well. 

The  application  of  swept  frequency  illumination  (Figure  6.7)  smoothed 
the  monochromatic  lobes  and  formed  a  uniform,  broad  backscatter  lobe. 


104 


The  effect  of  increasing  the  bandwidth  was  to  strengthen  the  off-center 
lobes  found  at  the  higher  and  lower  frequencies  and  hence  "fill"  the  sphere 
return.  Unfortunately  at  sweeps  above  about  +  10%  the  frequency  response 
curve  of  the  transducers  has  decreased  such  that  only  a  small  contribution 
is  obtained  from  the  extreme  frequencies  of  the  sweep#  hence  the  full 
effect  of  frequency  averaging  could  not  be  investigated.  This  is  especially 

evident  in  the  array  target  data . 

6.3.3  Array 

Figure  6.8  shows  the  acoustic  images  of  the  8- element  array  at 
three  monochromatic  frequencies.  The  images  were  recorded  for  several 
gain  settings  to  show  the  various  amplitude  levels  of  the  multiple  lobes. 

As  the  gain  is  increased  the  lower  level  lobes  become  evident.  Because 
of  the  gain  level  sensitivity  and  the  large  number  of  lobes  it  is  difficult 
to  determine  the  difference  in  the  lobe  structure  at  different  frequencies. 
Figure  6.9  compares  three  monochromatic  images  at  approximately  the 
same  loop-gain  levels#  i.e.  the  receiver  gain  is  set  to  compensate  for 

the  decrease  in  transducer  gain. 

Figure  6.10  shows  the  effect  of  adding  a  swept  frequency  signal 

of  two  different  sweep  widths:  +  100  KHz  and  +  50  KHz.  As  in  Figure  6.8# 
the  lobe  levels  were  interrogated  by  varying  the  receiver  gain.  The  +  50 
KHz  swept  signal  was  inadequate  to  completely  smooth  the  reradiation 
pattern#  however  the  +  100  KHz  swept  signal  provided  good  smoothing  at 
nearly  all  intensity  levels.  It  is  evident  that  a  greater  sweep  width  is 
required  to  completely  "fill"  the  ima^e. 
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monochromatic  frequencies . 


(a)  f=  1745  KHz 


(c)  f  =  1500 


i  ^  to  o  8  Effect  f  receiver  gain  variations  on  images  of  o-element 
'  array  at  three  monochromatic  frequencies 
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Figure  6 . 1 


(a)  Frequency  =  1640  +  100  KHz 


G  =  11  db  G  =  14  db 

(b)  Frequency  =  1640  +.50  KHz 


0  Effect  of  swept  frequency  on  image  :>f  8-element  array 
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CHAPTER  7 


CONCLUSIONS 

The  frequency  dependence  of  a  surface  of  known  statistics  was 
measured  and  compared  with  the  theoretical  predicted  results  of  a  theory 
based  on  the  Kirchhoff-Huygen's  Principle.  The  results  agree  over  a 
limited  frequency  range,  but  in  general  it  was  concluded  that  this  method 
does  not  adequately  predict  the  frequency  dependence  of  backscatter.  In 
addition  an  imaging  experiment  was  performed  employing  swept- frequency 
illumination.  The  experiment  demonstrated  the  improved  target  definition 
obtained  by  the  addition  of  frequency  averaging.  It  was  concluded  that  a 
panchromatic  system  should  offer  significant  advantages  over  monochroma¬ 
tic  systems  in  applications  such  as  terrain  imaging,  altimetry,  and 
tracking . 

Measurements  have  been  presented  for  the  frequency  dependence 

of  the  average  backscatter  power  from  a  surface  meeting  the  conditions 

of  the  Kirchhoff  method.  These  results  are  complemented  by  measurements 

of  the  frequency  dependence  of  the  average  backscatter  power  from  an 

extremely  rough  surface  for  which  the  Kirchhoff  method  is  not  applicable, 

but  which  is  indicative  of  many  natural  terrain  conditions.  When  the 

frequency  dependence  is  expressed  in  terms  of  the  wavelength  slope  as 
a 

X  ,  the  coefficient  a  was  found  to  vary  between  0  and  +  3.7  for  the 
fundamental  surface  and  -2  to  0  for  the  very  rough  surface.  The  measure¬ 
ment  frequency  range  was  approximately  6:1 .  A  comparison  of  the  mea¬ 
surements  results  with  the  results  predicted  by  the  Kirchhoff  method  was 
made  and  agreement  was  found  for  some  frequency  regions.  The  measure¬ 
ment  results  indicate  a  transi  donal  region  which  cannot  be  explained  by 
the  Kirchhoff  method.  In  addition,  the  measurements  support  previous 
suggestions  that  the  effective  statistical  parameters  of  the  surface  a^e 
both  angle  and  frequency  dependeint,  however  the  nature  of  this  dependence 
was  not  determined. 

Equipment  was  designed  and  constructed  to  implement  transmission 
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and  reception  of  a  broad -spectrum  acoustic  signal.  This  equipment  incor¬ 
porates  sufficient  flexibility  to  facilitate  examination  of  typical  radar  sys¬ 
tem  problems,  and  to  determine  the  effect  of  broad -spectrum  illumination 
on  the  reradiation  patterns  of  complex  targets.  The  system  was  used  to 
measure  the  decrease  in  positional  fading  and  the  averaging  of  complex  re¬ 
radiation  patterns  possible  with  the  incorporation  of  swept  frequency  illumi¬ 
nation.  The  technique  employed  to  obtain  panchromatic  illumination  was  to 
transmit  frequency  samples  of  a  continuous  frequency  spectrum.  Basic 
backscatter  measurements  were  recorded  using  a  high  frequency  sampling 
rate,  however  to  obtain  panchromatic  images  a  sample  rate  of  only  four 
samples  per  sweep  was  employed.  However  the  sampling  was  random  and 
hence  the  full  effect  of  panchromatic  illumination  was  achieved.  The  use 
of  panchromatic  illumination  substantially  improved  the  image  definition  of 
complex  targets . 

7.1  Recommendation  for  Further  Work 

The  experiment  resulted  in  defining  problems  requiring  further  work. 
The  most  important  of  these  are  the  following: 

A.  The  initial  concept  of  measuring  backscattered  power  over  a 
continuous  range  of  frequencies  was  hampered  by  the  nature 

of  the  acoustic  transducers  available.  In  addition  the  frequency 
dependence  of  the  target  did  not  warrant  so  fine  a  frequency 
sampling.  It  is  obvious  that  much  effort  is  required  before 
acceptable  transducers  are  available  for  such  experiments, 
however  it  does  not  appear  mandatory  that  they  operate  broad¬ 
band.  A  multiple  discrete  frequency  experiment  could  provide 
adequate  definition  of  the  frequency  dependence  of  most  targets 
except  for  conditions  to  be  discussed  next. 

B.  The  measurements  effectively  covered  a  frequency  range  of 
about  6:1 .  This  range  was  inadequate  to  enable  an  extension 

of  the  theory  to  fit  the  results.  The  analysis  of  the  results  does 
confirm  the  close  connection  between  angle  dependence  and 
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frequency  dependence  of  the  backscattered  power,  and  each 
appears  to  indicate  that  the  effective-parameters  concept  is 
appropriate.  To  determine  the  validity  of  the  suspected  effective- 
parameter  approach  or  the  ''size  filtering"  concept,  the  frequency 
range  must  encompass  several  transition  regions  of  structure 
sizes.  This  may  well  require  deliberately  including  structure 
sizes,  such  as  the  small  sand  particles,  that  violate  the  condi¬ 
tions  of  the  theory.  This  is  due  to  the  need  of  including  the 
regions  for  which  the  scatterers  are  approximately  the  size  of 
the  incident  wavelength.  To  illustrate  this,  note  that  the  fre¬ 
quency  dependence  curves  obtained  from  the  sanded  surface  did 
not  appear  to  have  reached  a  transition  region  which  showed  the 
existence  of  the  undulating  surface  (fundamental  surface) .  The 
sand  particles  dominated  the  return  throughout  the  measured 
range.  To  determine  the  transition  regions  may  require  a  very 
close  frequency  sampling,  This  realization  was  the  basis  of  the 
equipment  design,  however  the  restricted  frequency  range  covered 
did  not  warrant  the  capability  of  the  instrument.  The  eventual 
availability  of  a  wider  range  of  bet*?r  transducers  may  again  justify 
the  initial  design. 

C  ,  The  examination  of  the  Kirchhoff  method  showed  that  additional 
information  about  the  surface  statistics  needs  to  be  included 
or  at  least  a  new  interpretation  is  required  if  an  accurate  pre¬ 
diction  of  frequency  dependence  is  to  be  obtained.  The  present 
situation,  which  in  essence  says  that  the  greater  the  variance 
of  surface  heights  the  more  dominant  the  small  lag  distance 
region  of  the  autocorrelation,  is  contrary  to  physical  reasoning. 
Regardless  of  the  correlation  distance  associated  with  the 
.  urface ,  the  region  near  the  origin  of  the  autocorrelation  is 
dominated  by  the  small  structure  on  the  surface.  Obviously 
the  small  structure  does  not  contribute  to  the  fact  that  the 
surface  has  a  large  variance,  unless  the  surface  is  extremely 
abnormal.  Therefore  the  theory  forces  the  contradiction  that  the 
c  reater  the  variance  (a  large  structure  descriptor)  the  more  im- 
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portant  the  small  structure  becomes  in  determining  the  scatter. 

It  is  little  wonder  that  authors  are  lead  to  contend  that  the  high 
frequency  components  of  the  surface  have  more  influence  than 
the  low  frequency  components.  It  is  suggested  that  future 
work  on  this  problem  follow  the  work  of  Fung  and  Moore  (1964), 
which  was  a  composite  autocorrelation  function,  but  also  include 
variance  terms  associated  with  the  appropriate  subranges  of 
structure  sizes . 

D.  The  use  of  swept-frequency  illumination  in  image  recording  of 
complex  object  targets  and  area-extensive  targets  has  tremen¬ 
dous  potential  both  for  military  and  civilian  applications .  This 
is  subject  matter  for  an  extensive  research  program.  The  acoustic 
test  facility  developed  through  this  dissertation  effort  offers  an 
ideal  means  of  exploring  this  research  area.  However,  as  has 
been  mentioned  repeatedly,  the  full  potential  of  the  system 
cannot  be  realized  until  improved  acoustic  transducers  are 
developed . 
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APPENDIX  A 


GENERAL  KIRCHHOFF  METHOD 


This  is  a  summary  review  of  a  general  scattering  theory  based  on 
the  physical  optics  method.  The  treatment  follows  that  of  Beckmann  and 
Spizzichino  (1963)/  and  is  for  scalar  waves.  The  objective  is  to  deter¬ 
mine  the  predicted  frequency  dependence  of  backscatter  for  rough  surfaces. 

The  physical  optics,  or  Kirchhoff,  approach  has  been  employed 
by  several  authors  to  determine  scattering  from  both  "rough"  and  "smooth" 
surfaces.  The  general  approach  will  be  outlined  and  the  backscatter 
expressions  will  be  obtained  for  two  surfaces  with  different  autocorrela¬ 
tion  functions;  Gaussian  and  exponential. 

A.  1  General  Development 

The  scattering  geometry  with  which  we  will  be  concerned  is  shown 
in  figures  1  and  2. 


2 
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Figure  2 . 


The  scattered  field  at  the  observation  point  P  (figure  2)  is  given 
by  the  Helmholtz  integral. 


£,<?>  - 


where 


<T  e 


R' 


yM]  <u 

e>~/ 


(1) 


By  placing  P  in  the  Fraunhofer  zone  of  diffraction,  R’  —  <»#  then 
we  may  rewrite 


where  R  is  the  distance  from  P  to  the  origin. 

o 

In  (1)  we  are  concerned  with  the  values  of  E  and  3E/8n  which 
represent  the  field  and  its  derivative  normal  to  S .  The  Kirchhoff  or  physical 
optics  method  consists  essentially  of  approximately  the  values  of  E  and 
8E/8n  on  S  and  evaluating  (1) . 

We  will  now  make  an  assumption  which  is  fundamental  to  this 
approach.  It  is  assumed  that  the  field  at  any  point  on  the  surface  can  be 
approximated  by  the  field  that  would  be  present  on  the  tangent  plane  at 
that  point. 

With  this  approximation  the  field  on  S  will  be 

Es  =  ( /  +R)E,  (2) 


where  R  is  the  Fresnel  reflection  coefficient  of  a  smooth  plane.  Equation 
(3)  follows  from  (2)  by  differentiating  the  incident  and  reflected  waves . 


E -E, 
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(vf-  *)  ■  fte )  ,  J  f?)e  '*'Yl «) 

s 

note: 

(4a)  z/,"7  -  4.*  (4b) 

Employing  the  notation  H  ,  E  =  incident  field,  H. ,  E,  =  trans- 

00  ii 

mitted  field  and  H2 ,  E2  =  reflected  field  and  the  boundary  conditions: 

-yl  x  ( £0  +  )  *=  -yL  *  £, 

yi  X  ( fi0  +  f-/z  )  =  -rLx/^ 

then: 


where  0Q  =  angle  of  incidence  of  incident  wave.  Assuming  =  jjl2 ,  then 
from  (5) 
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where  6, 


=  i,e*  re9ion  2  is  free  space 


^  - 


CoS 


&0  -  l/~Y  ^  s/n  z  &o 


CoS  &o  —  t/Y  iS  •*-  1  @o  ' 
and  similarily  from  (6) 


/?,  ~  00-^  Y* L  s/jg*# 

Yzcos  &a+yry* ~ 


(7) 


(8) 


Now  returning  to  figure  2  and  using  (1) ,  (2)  and  (3)  we  compute 
the  scattered  field  at  a  point  P  in  the  far  field. 
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e, (p)  g'j[Ra,- o -  a.  *h)U « 


-Or 
c/s 


from  the  geometry,  figure  3  ,  we  can  write  the  unit  vector,  n,  as 


(10) 


The  ares,  ds ,  can  be  written  as 


(ID 
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The  quantity  [zS+z'+iJ*  accounts  for  the  projection  of  ds 
into  the  xy  plane.  Another,  more  illustrative,  way  of  writing  n  and  ds 
is  as  follows: 


■n  = 

cC  CoSyf  dy  sm ^  1 

= 

CoSyd  Cob  °C  ^ fiann*  *C  -i 

els  - 

d* 

clij  —  “ 

so 

CoS  oc 

Cospf 

~yL  ds 

~  [ 

Zx  +  Zu  *  & 

Substituting  (12)  into  (9) 

t>c  COSy4  eo&*X. 


c/xc/f 


Ez(p)  -  K  f  ( Rv-p  )-( Z v  A*  Oj  +&4)eJ  c/xc/y 


k,  -  l  6,  gL x  -  U  c&&  & , 

kt  =■  l  cos  @3  S/n  0Z  0lx  +  k  «S/n  &3  s/n  <2y  +  Jc  cos  Q±  d f 

V  =  ^7.  -  l £ " f  stn  B,  -  cos  Bj  3/n  Bz)a.x 

—  ( s/n  ^r/n  Bz)f,t  +  f-  cos  B,  -coiBz)^  J 
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p  -  lt  +  =  k jfsm  9,  +  CoS  &3S,n  &2)  a,x  V-(sin&jSin&z)  0.^ 

•*-  f -  cos j  Ot-*-  cos  J  ~j 

( Rv  ~  p)  (Zx  Ax  ■+/-£  A*  J  ~  (ft  ~  f  Jd xS*n  &,-(&+■  t  )  Cos  Sin  &z 

—  ( R+  I )  £tn  &3  Sirt.  Z y 

-(R-*  )  CoS  <9,  -(R+l)  Cos  Ot 


j  (k'-Lj).* 


cLx  dtc 


(14) 


b  *  —( R+/)  Sin  O3  Sir* 

C  =  —( R- 1 )  e*os  &z  -  1 )  cos  G, 

Assume  vertical  polarization  where  R  =  +1  (directly  applicable  to  acoustic 


-  2  cos  03  S/n  0Z  Zx 

7 


r 

d*  cl< 


(15) 


Now  we  will  attempt  to  eliminate  the  terms  Zx  and  Zy  from  the  integration, 
(15) .  Consider  the  term 


ja-iJ-R 

dxdfr 
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where 


2.  Cx.y.)  --  ft*,  %  > 


= 


e  4 ' 


iu. 


J  *  *  v* 


and  ^  e 


J'i^ejr'S) 


:  t/ 


and  integrate  by  parts 


£  -  K 


j* 


J'ii(x-3)  ;**+j*H  , 
-  c  u2 


e 


jtx+J'ji  */•«* 


dxda 


Consider: 


cli/  - 


T^y  C 


y  ?**<*.$) 


,  < ji  .,'•‘"'3^ 

.  /U  <^V, 


from  the  calculus 


•f dfii  -  /  7^>;  *  (Jx 


1/  » 


^  (x )  cl  x 


j  i/g  ?*C  x. 


c)  ?<*,  9-  c/if 

c)* 


(/  - 


J* 


X 


123 


We  now  assume  that  the  value  of  Z(x/y)  along  the  boundary  of  the  area  ds 
is  equal  to  the  mean#  i.e.  zero.  This  assumption  says  that  the  effect  of 
roughness  at  the  "edges"  is  zero  and  is  reasonable  valid  if  the  area  is 
large  relative  to  the  illumination  area.  Using  this  assumption#  the  first 
term  on  the  right  in  (16)  reduces  to 


j_  Jin  <SrX  3/n  K  Y 


(17) 


Equation  (17)  is  referred  to  as  the  "edge  effect  term"  and  is  gen¬ 
erally  small  enough  to  be  neglected  in  the  subsequent  development. 

In  the  same  manner  that  we  eliminated  Z  .  Z  can  be  eliminated 

X  y 

and  as  before  a  small  "edge  effect"  term  is  found.  The  total  edge  effect 
is  the  sum  of  these  two.  Now  (15)  can  be  written  as 

j  JcXo 

£z  (p)  -  -  /  g  _  j  /  -  Cos  sm  J$r 

Znf?o 


JT  Y 


where 


-  Sm  03  Sm  &z  %  -  Cos  C 

J 

m  /  [cm  0,  -  Cos  sm  J 

=»  “  A  ( Si  *  Si”  ) 

=■  -  A  { cos  6,  v  cos  &z  ) 


axd<£ 


(18) 


so 
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for  backs  catter. 


=  % 


-  TT 


(19)  becomes  (R  =  +l): 


Proceeding  with  (21)  as  before 
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(22) 


for  backs  natter, 


*  0%  J  Oj  —  7T  (22)  becomes  (R  =  -1): 


(23) 


Regardless  of  the  selection  of  P  =  +1  or  -1  the  power  calculation 
will  be  identical.  Therefore  the  proceding  development  does  not  represent 
the  vector  character  of  the  electromagnetic  field,  but  is  an  appropriate 
representation  of  the  acoustic  waves  being  considered. 


A. 2  Random  Rough  Surface 

We  will  now  consider  the  proceeding  as  applied  to  a  random  rough 
surface.  From  (19)  we  can  compute  the  mean  power  using  the  Poynting 
vector 


<Pj-i  Y'(eX> 


where  Y  =  admittence  of  free  space.  The  brackets  represent  the  ensemble 

®  -k 

average  of  the  terms  enclosed.  The  quantity  is 


clx  dg  clx'd.^' 


(24) 
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where  v-r  *  (k-kjr  has  been  expressed  as 

*  v*  X  +  ^3  + 

and  z  -  J  (x, y) 

Recall  that  it  was  assumed  that  <  £(x,y)  >  =  0  .The  surface 
term  £  is  the  random  variable  of  concern.  Equation  (24)  can  now  be 
written  in  this  form  of  an  average  as  follows. 

/  - *ZK //I f  e  (c  ydxd§dxdly{ 21 


i  rrr 


We  define  the  distance  between  the  points  (x,y)  and  (x' ,y')  as  t 
and  write  (25)  in  polar  coordinates. 


<E^)  *F>K  e 


}rcLrcLtf  (26) 


r  o 


Disregarding  the  limits  on  t  for  now,  we  note  that  the  integration  over 
0  is  the  zero  order  Bessel  function,  so  (2  6)  becomes: 


<^E*)-ZrrFtKjjo{ 


j^s-s'K 


71  e  /  ’•'air  (27) 


Consider  the  term  ^  C  )  ,  This  average  of  the  function 

^  r  i  n 


j  </*($-?•) 


p/7, 5')  Jfct 5' 


which  is,  by  definition  (Davenport  and  Root,  p.  52),  the  joint  character¬ 
istic  function  *(jvz,  -Jvz)  .  The  term  p(f,n  is  the  joint 

probability  distribution  of  the  random  variables  £  and  Using  the 
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characteristic  function  notation/  (27)  becomes 


(28) 


A. 3  Normally  Distributed  Surface 

We  will  consider  the  normal  distribution  of  surface  heights  and 
develop  the  backscatter  power  expression  to  determine  the  predicted 
frequency  dependence. 

The  characteristic  function  associated  with  the  normal  distribution 
is: 

X  (v*  )  s  e*P  (-£  &  \ 1 )  <30 

2 

where  a  -  variance  of  surface  heights .  The  characteristic  function 
associated  with  the  two  dimensional  distribution  is: 


=  etpf-vfar'O-  C)J 


(31) 


where  C  -  autocorrelation  function/ 
T  =  correlation  distance . 


Assuming  C  of  a  general  form  e 
expansion  results  in  the  following. 


-r'/T 


a. 

and  expressing  (31)  as  a  series 


•m  •  ° 


Following  Beckmann's  notation,  we  introduce  the  term  g  = 
convenience . 

Substituting  (32)  into  (2  8)  leads  to 


(32) 


(v±<r) 


'/» 


for 
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After  integration  (34)  becomes 

* 


)=  -rrf^Ke 


-^y4m 


(35) 


This  result  is  valid  for  g  less  than  unity.  For  g»l  (very  rough  surface) 

we  must  use  saddle  point  integration  to  reduce  equation  (34).  This  oper- 

-  r»/r  K 

ation  is  facilitated  by  expanding  the  autocorrelation  function  C  =  e 
into  a  MaClaurin  series  and  truncating  the  series  at  two  terms,  so  that 
(34)  reduces  to 


(£<  ft  *>  =  *»Pf- 

3  ( 


(36) 


Summarizing  and  rewriting  (35)  and  (36)  to  show  the  role  of  the 
wavelength,  X,  we  obtain: 


Case  1:  g«l  (slightly  rough  surface) 


TtTZFzKl  (RoF)  <fcexp^-2rro'X  Ycos  Ot 


T  fsm  Z0t  -ZstnO'  ton  9Z  COS  9^+Sth  %)j 37) 


X  exp£-TT 

Case  2:  g»l  (very  rough  surface) 


J 


129 


ZK.T 


/?  4rrorz  ( cos  B,  +  cos  Qz  )‘ 


x  q Xp  f  -  TZfsm  7 ~2sm  B,smB?cos  &K-b  S/n  ZQz)j 

L  4  o*  (  cos  B,  ~+  Cos  &i) *  7 

for  the  case  of  backs  catter,  -  Qz  ;  03  =  rr  ,  (37)  and  (38) 

become  respectively 

§■« / 

*)  =  '^Tr  Z>  exp ,/C  /<£  /7  V  \A~  *6>, 1 


» 


exp£-47rzTZ\  Z&tr\ZG,J 

/ 


/ EoEi)-  T  exp f-JL  iajn,1 9,1 

'  *  4/  ArrRSreos+e,  L2cr'  J 


(38) 


(39) 


(40) 


A. 4  Exponentially  Distributed  Surface 


In  this  section  the  treatment  of  A.  3  will  be  repeated  except  that 
the  autocorrelation  function  will  be  assumed  exponential,  i.e. 

C  =  exp  (-/  y!/t) 

The  characteristic  function  associated  with  the  two  dimensional 
exponential  distribution  is: 

\(vi  .-*>  )=  expf-v^cr'-f  <  -  e  /t)J  « 

which  may  be  written  as  follows 
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e  z 

Again  incorporating  •2r‘v*cr  and  substituting  (42)  into  (2  8)  we  obtain 


or 


(43) 


{ Eze/}  -  2rrFzKe  f^U^e 


-  rr>  IvI/t  , 

n'dr 


(44) 


After  integration  (44)  becomes 

<£-z£/)-~  ZnFzKc  -/  s 


\ 

-%  l  -  ™ 


-yYl 


T 


2. 


•V/'  s  / 

This  result  is  valid  for  g  less  than  unity,  for  g»l  (very  rough  surface) 
we  must  again  use  saddle  point  integration  to  reduce  equation  (44) . 
Proceeding  as  in  A.  3  equation  (44)  reduces  to  (Beckmann  1965) 


Tr>  if  f'n  2  ~t 


2.1  &Z 


(45) 


fa'*  (t*,)'. 


(46) 


Summarizing  and  rewriting  (45)  and  (46)  for  the  case  of  backscatter 
to  show  the  role  of  the  wavelength  .  X ,  we  obtain 

Case  1:  g«l  (slightly  rough  surface) 


t  e  e  *  y  -  77  3 ^-^2 

Z  *  PoZX[Xz+  TaJ6rr^Q]3/z 

e*p ’ F  rrZcrZCos> 


(47) 
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Case  2:  g»l  (very  rough  surface) 


Q  co$ z  o, 

i— -  — — —  ^ 

8rrf?JlfcoS*&l+QsmieJ  * 


(48) 
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APPENDIX  B 


GENERAL  AUTOCORRELATION  FUNCTIONS 


This  is  a  continuation  of  the  discussions  presented  in  section  A. 3 
and  A. 4  of  Appendix  A  with  an  elaboraf  on  relevent  to  the  comparison  of 
measured  and  predicted  results  presented  in  section  5.4. 

In  section  4.3  it  was  shown  that  the  autocorrelation  of  the  funda¬ 
mental  surface  was  fit  by  the  Gaus:,ian  function  for  the  region  away  from 
the  origin,  and  by  the  following  function  for  the  region  near  che  origin. 


C(r)  =  exp 


(B.l) 


Substitution  of  this  autocorrelation  function  into  the  Kirchhoff  ex¬ 
pression  (equation  2.10)  poses  considerable  difficulty  when  attempting  the 
integration,  therefore  the  following  approach  is  employed  (Fung  and  Leovaris, 
1968).  Recall  that  to  obtain  (A. 39)  and  (A. 40)  the  autocorrelation  function 
C (t)  was  expanded  about  t  =  0  and  the  series  was  approximated  by  the  first 
two  terms.  The  first  derivative  C'(t)  is  necessarily  zero  since  C(t)  is 
analytic,  hence  the  autocorrelation  function  is  approximated  by 


C  (r) 


/  *  c  To)  r  * 

Z 


Consider  that  instead  of  expanding  about  t  =  0,  the  expansion  is  made  about 
t  =  T-^ ,  where  is  the  value  for  which  the  integrand  (equation  2.10)  .s  a 
maximum.  Then  for  small  t  the  approximation  is 

CM  -  /  +  C  ?rt)  y 
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* 


where  C'(tj)  is  the  slope  of  the  autocorrelation  function  at  . 

Using  this  approach  the  equivalent  expression  to  (A. 40)  is 
(omitting  a  constant): 


(B.2) 


r 


Equation  (B.2)  is  identical  in  form  to  equation  (A. 48)  obtained  for 
the  exponential  autocorrelation.  However,  using  (B.2)  allows  the  use  of 
many  autocorrelation  functions  for  which  equation  (2.10)  cannot  other¬ 
wise  be  readily  integrated. 

The  first  derivative  of  the  autocorrelation  function  in  equation 


(B.3) 


If  is  selected  to  be  a  constant,  the  frequency  and  angular 
dependence  of  (B.2)  is  identical  with  that  of  the  exponential  autocorrelation 
shown  in  section  5.4.  Letting  be  the  point  for  which  Texp  [  -g  (1  -  C(t))  ] 
in  equation  (2.10)  is  maximum  leads  to  values  for  equation  (B.3)  as  listed 
in  Table  B.l,  using  the  profile  measurements  of  the  fundamental  surface. 
Figure  B.l  shows  <j°  versus  wavelength  using  equation  (B.2)  and  the 
values  of  C'^)  shown  in  Table  B.l . 
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Using  this  approach,  Fu  -g  and  Leovaris  (1968)  determined  values 
of  Q  which  provided  a  reasonable  fit  to  experiment  results  of  angle  de¬ 
pendence  at  five  wavelengths.  From  these  values  of  Q,  the  slope  term, 
C'(t^)  was  calculated  and  was  found  to  compare  favorably  to  the  measured 
slope  of  the  surface  autocorrelation.  Their  work  shows  that  the  variation 
of  angular  dependence  with  wavelength  is  reasonably  well  predicted  by 
equation  B.2.  This  method  is  analagous  to  employing  an  effective  para- 
meter,  T/cr  ,  which  varies  with  wavelength.  Using  the  expression  assumed 
in  equation  5 . 3 

f?e  Rf(G)^- A j  (5.3) 

2 

where  R  =  effective  parameter,  (T/ct  ) 

C 

2 

R  =  actual  parameter,  (T /a  ) 


then  from  Table  B.l  the  value  of  j  is  approximately  -0.5.  This  results 
in  a  predicted  a  for  0°  incidence  of  1 .0  which  compares  with  the  measured 
o?=  1.03  (below  1.25  MHz)  and  a=  2.78  (above  1.25  MHz).  The  predicted 
a  for  15°  and  30°  incidence  is  also  in  good  agreement  with  the  measure¬ 
ments  of  a  below  1.25  MHz.  However,  the  data  above  1.25  MHz  is  not 
fit  by  the  theory. 
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TABLE  B.l 


X  (mm) 

©  (degrees) 

* 

(mm) 

C,(t1) 

0.5 

0 

0.508 

-.0080 

15 

0.508 

-.0080 

30 

0.508 

-.0080 

45 

1 .016 

-.0113 

60 

1.016 

-.0113 

1.0 

0 

1.016 

-.0113 

15 

1.016 

-.0113 

30 

1.524 

-.0138 

45 

1.524 

-.0138 

60 

2.540 

-.0175 

1.5 

0 

1.524 

-.0138 

15 

1.524 

-.0138 

30 

2.032 

-.0158 

45 

2.540 

-.0175 

60 

3.556 

-.0203 

3.0 

0 

3.556 

-.0203 

15 

3.556 

-.0203 

30 

4.064 

-.0215 

45 

5.588 

-.0243 

60 

** 

— 

4.74 

0 

6.604 

-.0256 

15 

7.112 

-.0262 

30 

9.655 

-.0281 

45 

** 

— 

60 

** 

— 

These  values  correspond  to  the  nearest  point  for  which  a  measurement 
was  recorded. 

The  calculated  value  of  exceeded  a  meaningful  value. 
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Figure  B .  1 

PREDICTED  FREQUENCY  DEPENDENCE 
Varying  Slope  Method 
Fundamental  Surface. 
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